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This  paper  presents  a  new  methodology  for  assessing  the  performance  of 
area  defense  weapon  systems  against  airbreathing  strategic  threats.  The 
methodology  is  probabilistic  rather  than  simulative.  It  intertwines  sur¬ 
veillance  sensor,  interceptor,  and  command,  control  and  communications  net¬ 
work.  Furthermore,  it  allows  for  incorporating  the  propensity  of  the  air- 
breathing  threats  to  maneuver.  The  probabilistic  approach  provides  the  tools 
for  direct  sensitivity  analyses  against  overall  weapon  system  performance 
including  tracking  errors.  The  approach  is  amenable  to  direct  interceptor 
scheduling  and  yields  closed-form  solutions,  thus  permitting  relatively  easy 
(and  inexpensive!)  programming. 

The  analytical  presentation  of  the  model  is  supplemented  by  a  computeriza 
tion  and  a  numerical  example. 
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the  manuscript,  and  Mr  Jack  Bennett  for  drawing  the  figures. 

A  J-ii  ■*£  / 


COURTLANDT  L.  BOHN,  ILt,  USAF 
Staff  Physicist 

Directorate  of  Aerospace  Studies 


This  report  has  been  reviewed  and  is  approved  for  publication. 


BRUNO  J.  MANZ,  Assistant  Chief 
Studies  &  Evaluation  Division 
Directorate  of  Aerospace  Studies 


PAUL  J./VALLERIE,  Lt  Col,  USAF 
Director  of  Aerospace  Studies 
DCS/Plans  and  Programs,  HQ  AFSC 


This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is 
releaseable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 


Accession 

KTIS  GP.Ail 
DTIC  TAB 
U*’anrouiiced 
Justification — 


9 


By. 

Di 


I 


t.rieutlon/ 

■  i.'.-'Mlity  Code 
and /or 
Upecial 


r 


Diit 


UNCLASSIFIED _ 

StCUMITV  CLASSIFICATION  OF  THIS  PASS  f'Whan  Omtm  Enttrtd) 


REPORT  DOCUMENTATION  PAGE 


lv  REPORT  NUMBER 


BAfP  READ  instructions 

rAV,c _ BEFORE  COMPLETING  FORM 

T  GOVT  ACCESSION  NO.  1-  RECIPIENT'S  CATALOG  NUMBER 


«•  title  Im dSuttttimy  - - 

/V  AIRBORNE  PLATFORMS  VERSUS JIRBREATHING 

!j  )*% TRATEGIC">JHREATS./'A  performance  assessment 

^ '  WfTHODOLOGT,  ^  *  ' 


/'?]  Final  h}'tpTi 
I  fJ\  Jan  -  Jun  498 Si 


UOC.COVEPED 


Tr-AuTH«!V.; . -  - . I  ■ 

^tet-Courtlandt  L/Bohn 

sRFA.  Foster/Cooper 

-^e  Bruno  JyManz  _ 

Sr  PERFOHBISTG^Oire ANIMATION  NAME  ANO-AOORESS 

Directorate  of  Aerospace  Studies  / 
Kirtland  AFB,  NM  87117  7 


I')  1 


••  CONTRACT  or  GRANT  NUMBCRf*; 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


1 1  1.  CONTROLLING  OFFICE  NAME  AND  AOORESS 


^  U.  RE-PORT  PATE 

//,  June  198j6  / 


IT  MONITORING  AGENCY  NAME  A  ADDRESS!"//  dllfmrmnt  Iron  Controlling  Of/ieoJ  15.  SECURITY  CLASS,  (ot  thi,  rmpon) 

UNCLASSIFIED _ 

ts«.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

<6.  DISTRIBUTION  STATEMENT  (o I  thi,  Rmpon) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetrect  entered  In  Block  20,  It  different  from  Report) 


t».  SUPPLEMENTARY  NOTES 


I  VS.  KEY  WORDS  (Continue  on  reveree  tide  if  neceeeery  end  Identify  by  block  number) 


Strategic  Defense 
Airbreathing  Threats 
Airborne  Platforms 


Applied  Probability  Theory 
Alr-to-Air  Engagement 
Distribution  Functions 
Homing  Simulation 


[ABSTRACT  ( Continue  on  reeeree  a Ida  If  neceeeary  end  Identify  by  block  number) 

(This  paper  presents  a  new  methodology  for  assessing  the  performance  of  area 
defense  weapon  systems  against  airbreathing  strategic  threats.  The  method¬ 
ology  is  probabilistic  rather  than  simulative.  It  intertwines  surveillance 
sensor.  Interceptor,  and  command,  control  and  communication  network. 
Furthermore,  It  allows  for  incorporating  the  propensity  of  the  airbreathing 
threats  to  maneuver.  The  probabilistic  approach  provides  the  tools  for 
direct  sensitivity  analyses  against  overall  weapon  system  performance 
including  tracking  errors.  The  approach  is  amenable  to  direct  Interceptor 


DD 


EDITION  OF  I  NOV  «S  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  fiflwn  Ow« 


=3 

fO  h%  '  '/ 


CUWITV  CLASSIFICATION  OF  this  PAOZOTh M  Omit  Sntmrmd) 


scheduling  and  yields  closed-form  solutions,  thus  permitting  relatively  easy 
(and  inexpensive)  programming. 

The  analytical  presentation  of  the  model  is  supplemented  by  a  computerization 
and  a  numerical  example. 


UNCLASSIFIED 


ItCuaiTY  CLASSIFICATION  OF  THIS  F*GE'Wh«n  Dim  En r#r*tf> 


MODEL  OVERVIEW 


This  model  is  designed  to  support  an  analysis,  the  objective  of  which  is 
to  determine  the  potential  of  candidate  air  defense  engagement  systems  to 
counter  a  Soviet  airbreathing  attack  against  the  continental  United  States  in 
the  1985-2010  time  frame.  One  class  of  engagement  systems  centers  on  area 
defense  with  an  air-based  weapon  platform.  This  is  the  class  the  Airborne 
Platforms  Versus  Airbreathing  Strategic  Threats  (APVAST)  model  treats.  A 
possible  measure  of  merit  centers  on  expected  number  of  penetrators  negated, 

<  Nk>  .  By  combining  salient  features  of  interceptor,  sensor,  penetrator, 
and  command,  control,  and  communications  (C^),  the  APVAST  methodology  leads 
to  a  calculation  of  expected  Nil!  probability,  <PK>  ,  from  the  conventional 
kill  probability,  P^,  which  in  turn  gives  <  N^>  .  The  methodology  is 
probabilistic  rather  than  simulative  and  is  adaptable  to  either  an  abstract 
treatment  of  scenarios  or  to  a  full-blown,  concrete  campaign  model. 

A  standard  approach  to  modelling  air  defense  problems  involves  Monte  Carlo 
techniques.  It  requires  as  the  initial  input  a  concrete  campaign  laydown, 
wherein  individual  penetrators  are  explicitly  assigned  to  individual  targets. 
By  exercising  the  model  many  times,  the  user  can  calculate  expected  values  and 
corresponding  distribution  functions  directly.  While  the  Monte  Carlo  tech¬ 
niques  are  straightforward  and  easy  to  understand,  they  invariably  lead  to 
inefficient,  i.e.,  time-consuming,  software.  In  turn,  sensitivity  analyses 
become  expensive. 

These  drawbacks  are  altogether  avoided  using  probabilistic  techniques. 
These  techniques  are  analytically  demanding  and  less  intuitive,  but  they  lead 
to  a  model  which,  when  exercised  only  once,  is  equivalent  to  an  infinite 
number  of  Monte  Carlo  runs.  The  model  can  accept  abstract  weapons  laydown  as 
an  input,  and  it  also  can  accept  a  concrete  laydown.  It  is  a  tool  for  doing 
immediate  sensitivity  analyses. 

The  basic  APVAST  goal  is  to  convert  conventional  kill  probability  into 
expected  kill  probability.  The  principle  for  doing  so  is  elementary.  Every 
sensor  reading  is  subject  to  uncertainties.  Since  a  sensor  measures  both 
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penetrator  position  and  velocity,  there  are  six  uncertainties,  one  for  each  of 
the  three  position  degrees-of-freedom  and  one  for  each  of  the  three  velocity 
degrees-of-freedom.  With  each  of  the  six  uncertainties  is  associated  a 
corresponding  distribution  function.  For  example,  if  the  uncertainties  are 
all  Gaussian,  then  each  of  the  associated  distributions  is  Gaussian,  and  the 
penetrator  is  "centered"  within  a  six-dimensional  Gaussian  ellipsoid.  (The 
ellipsoid  set  up  by  the  three  position  uncertainties  is  called  an  "error 
basket.")  As  long  as  a  penetrator  remains  within  sensor  coverage,  the  ellip¬ 
soidal  axes  are  determined  solely  from  sensor  parameters  depending  on  the 
penetrator  spacing.  However,  once  the  penetrator  leaves  sensor  purview,  its 
position  uncertainty  grows  with  time  and,  consequently,  the  corresponding 
expected  kill  probability  shrinks  through  the  growing  uncertainty  inherent  in 
the  interceptor's  vectoring  information. 

Two  other  factors  affect  the  accuracy  of  the  vectoring  information.  One 
3  3 

is  the  C  network.  For  example,  the  C  network  may  allow  for  continuous 
update,  periodic  update,  or  no  update  of  the  interceptor's  vectoring  instruc¬ 
tions,  and  the  <  PK>  shrinks  correspondingly.  As  a  case  in  point,  if  the 
C3  cannot  allow  for  vectoring  update,  then  intercept  is  driven  solely  by  the 
last  available  sensor  reading.  The  corresponding  error  basket  will  have  grown 
from  the  time  of  the  last  sensor  reading  up  to  the  attack  time  with  the  con¬ 
comitant  shrinking  of  <  PK>  ,  just  as  it  would  when  the  penetrator  leaves 
sensor  coverage.  Another  important  C3  property  is  inherent  dead  time.  For 
non- zero  dead  time,  the  vectoring  instructions  are  based  on  "dated"  sensor 
readings. 

The  other  important  factor  affecting  the  error  basket  is  penetrator 
maneuvers.  Once  sensor  readings  become  unavailable  and  the  error  basket 
begins  to  grow,  that  growth  is  further  exaggerated  by  penetrator  maneuvers. 

The  error  basket  becomes  "warped",  i.e.,  no  longer  Gaussian,  and  because  it  is 
larger,  penetrator  maneuvers  serve  to  further  degrade  <  PK>  .  Though  the 
calculation  is  difficult,  the  warped  error  basket  can  be  analytically  deter¬ 
mined  through  a  series  of  convolutions  (integrations)  over  the  sensor- induced 
and  maneuver- induced  distribution  functions. 
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The  <Pj<>  calculation  against  a  given  penetrator  also  hinges  on  deter¬ 
mining  the  expected  attack  time,  <  tA>  ,  and  corresponding  expected  attack 
position,  <  rA>  ,  to  which  the  interceptor  flies  to  attack  the  penetrator 
based  on  its  vectoring  instructions.  This  calculation  impacts  PK  which,  in 
turn,  impacts  <  PK>  .  If  the  penetrator  does  not  maneuver,  then  <  t^>  , 

<rA>  arise  through  elementary  vector  analysis.  If,  on  the  other  hand, 
the  penetrator  does  maneuver,  then  (given  vectoring  information  update)  these 
maneuvers  drive  the  interceptor  to  maneuver  and  home  for  kill.  The  proba¬ 
bilistic  approach  offers  a  natural  framework  within  which  to  account  for  this 
homing. 

In  summary,  an  alternative  to  Monte  Carlo  techniques  is  now  available  for 

air  defense  modelling.  The  key  is  a  probabilistic  approach  through  which 

3 

interceptor,  sensor,  C  ,  and  penetrator  properties  can  be  immediately  tied 
together  to  provide  a  completely  analytical  solution  for  <PK>  •  At  user's 
discretion,  the  methodology  can  be  applied  either  abstractly  or  in  a  concrete 
campaign  model,  and  its  flexibility  allows  for  sensitivity  analyses  to  both 
performance  parameters  and  scenarios. 
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I.  INTRODUCTION 


The  defense  of  our  land-based  strategic  forces  against  the  airbreathing 
leg  of  the  Soviet  triad  is  a  major  driver  of  current  Air  Force  development 
planning  activities.  One  manifestation  is  the  November  1979  tasking  of  the 
Directorate  of  Aerospace  Studies  (office  symbol  AFCMD/SA)  by  Headquarters  Air 
Force  Systems  Cotrmand  (AFSC)  to  support  an  AFSC  Air  Defense  Mission  Analysis. 
The  major  analysis  objective  is  to  determine  the  potential  of  candidate 
surveillance  and  engagement  systems  to  protect  the  continental  United  States 
from  the  1985-2010  Soviet  airbreathing  threat.  For  the  usual  cost- 
effectiveness  assessment,  it  is  first  necessary  to  quantify  the  effectiveness 
of  a  large  number  of  systems  and  system  combinations.  For  the  initial 
screening,  AFCMD/SA  has  chosen  as  measure  of  merit  the  number  of  penetrators 
killed, <  Nk>.  The  objective  is  all-encompassing;  sensors  could  be  in  any 
basing  mode  from  ground-based  through  space-based  and  could  rely  on  any  one  of 
several  advanced  technologies,  while  far-term  interceptors  could  be  designed 
for  either  area  defense  or  terminal  defense,  or  could  have  a  layered  defense 
profile,  and  could  also  be  in  any  basing  mode.  This  paper  concentrates  on 
area  defense  systems,  one  leg  of  which  is  constituted  by  mobile  airborne 
weapon  platforms.  It  presents  a  methodology  for  quantifying  the  performance 
of  the  entire  defense  system  consisting  of  sensors,  interceptors,  and  command, 
control,  and  communications  (C  )  capabilities. 

The  generic  scenario  under  study  is  represented  in  figure  l._  A  raid  of 
strategic  airbreathing  penetrators  originates  from  home  bases,  penetrates  a 
tactical  warning  fence,  proceeds  to  the  targets,  and  then  heads  for  recovery 
bases.  Their  goal  is  to  attack  all  targets  at  the  same  time  to  avoid  early 
tactical  warning  induced  by  early  weapon  detonation.  Each  raid  member 
penetrates  from  a  highly  favored  direction;  for  example,  the  most  favored 
direction  of  a  bomber  raid  penetrating  the  Distant  Early  Warning  Line  is 
obviously  from  north  to  south.  The  distance  from  the  home  base  to  the  point 
of  fence  penetration  is  denoted  Rgpj  the  distance  from  the  point  of  fence 
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Home/Staging 


Figure  1.  Schematic  Outline  of  the  Problem 
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penetration  to  the  target  is  RpT;  and  the  distance  from  the  target  to  the 
recovery  base  is  R^.  Clearly,  if  the  individual  penetrator  is  to  recover 
safely,  then  the  sum  of  Rgp,  Rpy,  and  R-j-^  cannot  exceed  the  total  range 
capability  Rp  of  the  penetrator.  Impeding  the  attack  is  an  active  defense 
system  consisting  of  the  tracking  sensor,  C  network,  and  interceptor. 
Generally  these  elements  are  separate;  however,  in  certain  revolutionary 
concepts  they  may  be  collocated,  in  which  case  the  interceptor  becomes 
autonomous.  This  "trio"  of  subsystems  comprising  the  total  defense  system 
presents  a  rather  difficult  modelling  problem  which  is  further  complicated  by 
the  fact  that  an  airbreathing  penetrator  can  maneuver  by  changing  ’ts  heading, 
or  its  speed,  or  both. 


A  Monte  Carlo  approach  would  lead  to  an  excessive  number  of  cases,  while 
the  simple  simulation  of  one  "representative"  case  would  be  so  strongly  biased 
that  the  results  would  be  worthless.  Both  pitfalls  are  avoided  by  the 
probabilistic  approach  which,  however,  is  analytically  demanding. 


The  essence  of  the  probabilistic  approach  is  brought  out  by  the  following 
considerations.  A  sensor  reading  indicating  position,  speed,  and  heading  has 
associated  errors,  or  standard  deviations.  As  the  penetrator  proceeds  to  the 
target,  these  errors  cause  the  position  uncertainty  to  grow  in  a  quantifiable 
way;  maneuvers  worsen  the  position  uncertainty.  The  magnitude  of  the  position 
error  at  the  intended  time  of  intercept  also  depends  on  C  capabilities, 
that  is,  on  the  dead  time  between  sensor  reading  and  relaying  of  vectoring 
information.  Finally,  the  expected  probability  of  kill,  <PK>,  depends  on 
both  the  position  uncertainty  at  the  intended  time  of  intercept  and  on  the 
interceptor  capabilities.  Given  fixed  interceptor  capabilities,  a  larger 
position  uncertainty  results  in  a  smaller  <P^>. 

For  optimal  scheduling  of  the  mobile  intercept  force,  it  is  necessary  to 
know  the  times  at  which  and  the  positions  from  which  the  penetrators  are 
attacked.  Moreover,  to  insure  that  a  penetrator  is  killed  before  it  can 
attack  its  assigned  target,  one  also  needs  to  know  the  time  and  position  of 
kill.  These  requirements  are  met  by  an  approach  which  models  interceptor 
homing.  Of  course,  this  approach  has  to  account  for  penetrator  maneuvers. 

All  of  these  requirements  are  satisfied  in  the  most  natural  way  by  the 
probabilistic  approach. 


The  following  sections  develop  the  probabilistic  methodology.  Section  II 

develops  the  necessary  tools  for  quantifying  the  growth  of  both  sensor-induced 

and  maneuver- induced  position  uncertainties,  and  it  discusses  how  to  use  these 

tools  to  calculate  <PK>.  Section  III  lays  out  the  basic  equations  which 

couple  attack  geometry  and  interceptor  capabilities  without  allowing  for 

maneuvers.  The  tools  developed  in  section  II  and  the  basic  equations  of 

section  III  form  the  foundation  for  section  IV  which  folds  in  penetrator 
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maneuvers  to  describe  interceptor  homing,  subject  to  C  capabilities. 

Section  V  provides  an  example  to  demonstrate  the  model's  adaptability;  it 
describes  modelling  of  an  air  defense  system  which  centers  on  standoff  guided 
air-to-air  missiles  (AAMs)  --  an  at  least  conceptually  possible  alternative. 
The  theory  concludes  with  section  VI  which  sketches  the  logical  steps  required 
to  determine  expected  number  of  penetrators  negated.  Tedious  calculations  are 
relegated  to  the  appendices.  Appendix  A  discusses  the  transformation  of 
standard  deviations  from  one  coordinate  system  to  another.  Appendix  B 
presents  an  approximate  solution  of  an  integral  surfacing  in  section  II. 
Appendix  C  presents  the  APVAST  program  and  serves  as  an  abbreviated  user's 
manual.  The  reader  who  needs  only  to  know  how  to  the  use  the  program  is 
referred  directly  to  Appendix  C. 

The  model  has  been  computerized.  Appendix  C  offers  the  listing  of  the 
computer  program,  a  list  of  translation  of  the  symbols  of  the  theory  into 
computer  language,  and  a  numerical  example  in  input-output  form. 
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II.  PROPAGATION  OF  POSITION  UNCERTAINTY 


Measurements  are  always  subject  to  uncertainty  and  are  generally  associa¬ 
ted  with  distribution  functions.  In  the  present  problem,  penetrator  position 
uncertainty  arises  initially  due  to  sensor  imperfections.  As  time  passes  on, 
ana  the  penetrator  proceeds  to  the  target,  the  position  uncertainty  evolves 
due  to  three  factors  —  penetrator  motion,  sensor  motion,  and  penetrator 
maneuvers.  This  section  develops  the  tools  needed  to  quantify  this  evolu¬ 
tion.  It  begins  with  a  discussion  of  the  distribution  functions  pertaining  to 
the  initial  measurements.  It  then  calculates  the  evolution  of  the  final  posi¬ 
tion  distribution  functions.  Thereafter,  it  turns  to  the  use  of  the  position 
distribution  functions  in  calculating  <P«>- 

A.  INHERENT  DISTRIBUTION  FUNCTIONS 

i.  Problem  Schematics.  The  geometry  of  the  air  defense  problem  pro¬ 
vides  a  natural  framework  for  characterizing  the  inherent  distribution  func¬ 
tions.  Figure  2  shows  an  overhead  view  of  the  problem  geometry  assuming  a 
fiat  earth.  The  justification  for  neglecting  the  earth's  curvature  lies  in 
the  fact  that  the  engagement  of  penetrators  by  sensors  and  interceptors,  due 
to  the  limited  range  cf  the  latter,  always  takes  place  within  dimensions  which 
are  small  compared  to  the  earth's  radius.  In  figure  2,  an  interceptor  air 
base  CAB j  serves  as  the  origin  for  the  cartesian  reference  system  (x,y,z). 

The  y-axis  is  aligned  toward  the  most  favored  direction  from  which  the  raid  is 
coming.  For  concreteness,  that  direction  is  taken  to  be  north.  Centered  on 
the  sensor  (S)  is  the  sensor-based  cartesian  system  (x$»Ys,zS^' 

Also  shewn  in  figure  2  is  the  tactical  warning  fence  F.  It  separates 
the  defended  from  the  undefended  airspace.  The  basic  assumption  is  that  the 
penetrator  will  not  maneuver  before  it  reaches  the  fence,  but  it  generally 
will  maneuver  thereafter.  The  penetrator  trajectory  originates  from  the  home 
base  PH>  penetrates  the  fence  at  the  point  Pp,  and  proceeds  to  the  target 
P.p  Whereas  the  trajectory  section  from  P^  to  Pp  is  straight  by  assump¬ 
tion,  the  section  from  Pp  to  PT  is  generally  curved  due  to  maneuvers.  In 
figure  2,  the  curved  path  from  Pp  to  PT  is  an  example  of  an  actual  pene¬ 
trator  trajectory. 
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At  some  time  t$.  the  sensor  peforms  an  observation,  or  measurement. 

At  this  time,  the  penetrator  occupies  some  postion  P ( t^ )  on  the  curved  path, 
and  has  some  velocity  v(tc) ,  as  is  shown  in  figure  2.  The  sensor-determined 
penetrator  position  and  velocity  are  P(t^)  and  vU^),  also  shown  in 
figure  2.  Generally,  due  to  sensor  imperfections,  P ( )  and  v(t^)  will 
deviate  from  P(t^)  and  v(t^). 

It  is  now  important  to  understand  the  following.  If  somebody  measures 
an  unknown  quantity,  say  x,  and  quotes  the  result  of  his  measurement,  he 
always  quotes  that  value  of  x  which,  in  his  mind  and  with  regard  to  his 
information,  is  the  "most  probable"  value  of  x.  If  we  denote  this  value  by 
<x|Inf>,  then  we  may  say  that  this  is  the  result  of  the  measurement  which 
yielded  the  information  listed  in  the  condition  compartment. 

A 

These  conditions  apply  to  the  sensor-observed  data  P(t^)  and 
v(t<j).  They  are,  from  the  point  of  view  of  the  sensor,  the  most  probable 
penetrator  position  and  velocity  at  the  time  ts. 

2.  Perturbation  Approach.  On  the  trajectory  section  from  Pp  to 
Py,  the  penetrator  generally  maneuvers.  Therefore,  the  trajectory  section 
is  generally  not  straight.  However,  in  the  absence  of  any  information  about 
the  actual  maneuvering,  the  expected,  or  average,  trajectory  is  a  straight 
line.  This  expected  trajectory  is  shown  in  figure  3,  where  it  is  labeled  L. 
Clearly,  L  is  the  continuation  of  the  straight  line  from  PH  to  Pp. 

The  gist  of  the  perturbation  approach  is  the  following.  The 

A 

sensor-observed  position  P(t^)  is  the  most  probable  penetrator  position  from 
the  point  of  view  of  the  sensor,  but  not  from  the  point  of  view  of  the 

A 

analyst.  In  fact,  from  the  analyst's  point  of  view,  P(t^)  is  simply  an 
unknown  sample  value.  What,  then,  is  the  most  probable  sensor-determined 
penetrator  position  from  the  point  of  view  of  the  analyst?  Here  the  straight 
line  L  from  figure  3  comes  into  play.  This  is  the  expected,  or  average, 
trajectory.  From  the  analyst's  point  of  view,  it  is  the  most  probable 
trajectory.  Therefore,  from  the  analyst's  point  of  view,  the  most  probable 
sensor-determined  penetrator  position  at  some  given  time  tQ  lies  on  L. 
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It  is  for  these  reasons  that  the  perturbation  approach  replaces  the 

A 

unknown  sample  value  of  the  sensor-determined  penetrator  position,  P(t<J,  by 

A  w 

a  suitably  selected  point  <P(tQ)>  on  L.  This  is,  from  the  analyst's  point 
of  view,  the  most  probable  sensor-determined  penetrator  position  at  the  time 

V 

From  now  on,  the  "hat  symbol"  (~)  will  be  used  to  indicate  sensor- 
determined  or  sensor-dependent  values,  as  distinct  from  maneuver- induced 
values  to  be  introduced  later.  The  ^-values  are  conditional.  The  condition 
is  "no  maneuver."  The  "bracket  symbol"  (<  >)  will  always  carry  the  meaning 
"most  probable  from  the  analyst's  point  of  view." 

In  principle,  the  penetrator  has  six  degrees-of-freedom  —  three  in 
position  and  three  in  velocity.  At  the  outset,  however,  we  assume  that  the 
penetrator  always  flies  within  a  horizontal  plane,  thus  eliminating  one 
velocity  degree-of-freedom.  This  assumption  is  generally  valid  except  near 
the  target.  Each  penetrator  now  has  five  degrees-of-freedom,  where  the 
velocity  is  now  represented  by  the  penetrator  speed  v  and  heading  <j>. 

(Heading  will  always  be  taken  from  north,  as  in  figure  3.) 

The  point  <P(tQ)>  on  the  expected  path  L  must  be  "suitably 
selected."  In  fact,  the  coordinates  of  <P(tQ) >,  i.e.,  [<x(tQ)>, 

<y ( t Q ) > ] »  are  both  functions  of  past  maneuvers  insofar  as  these  maneuvers 
reduce  the  "effective  speed"  of  the  penetrator  along  its  expected  path.  In 
essence,  then,  we  model  maneuver  capability  in  part  through  the  calculation 
of  <P(t  )>.  Better,  we  model  maneuver  capability  by  treating  it  as  a 
perturbation  of  the  penetrator' s  movement  along  its  expected  path.  This  is 
the  essence  of  the  "perturbation  approach." 

A 

Connected  to  the  most  probable  observed  penetrator  position  <P(trt)> 
is  the  most  probable  observed  velocity  <v>  which  is  shown  in  figure  3.  We 
decompose  this  velocity  into  speed  and  heading.  The  most  probable  observed 
speed  is  <v>  which  is  just  the  operational,  or  cruising,  speed  of  the 
penetrator.  It  is  a  fixed  parameter;  that  is,  it  keeps  the  same  value  at  all 
times.  The  most  probable  observed  heading  <(£>  is  just  the  heading 
characterizing  the  expected  path  L.  Because  L  is  straight,  <ij>  assumes  the 


same  value  along  all  points  of  L.  Finally,  we  choose  as  the  most  probable 
observed  penetrator  altitude  <z>,  the  operational,  or  cruising,  penetrator 
altitude  h. 


The  geometry  of  the  perturbation  approach  is  illustrated  in  figure  4, 
which  focuses  on  that  section  of  the  expected  penetrator  trajectory  L  which 

A 

contains  the  most  probable  position  <P(tn)>.  The  penetrator's  freedom  to 
maneuver  |s  now  accounted  for  by  the  "change-of -speed"  v  and  the  "change  of 
heading"  iji.  From  now  on,  the  "plus  symbol"  (+)  will  be  used  to  identify 
maneuver- induced  or  maneuver- dependent  values. 

It  now  follows  that  the  "state  of  the  penetrator"  at  the  time  t  is 

A  A  A  a  0 

described  by  seven  distribution  functions.  These  are  X[x(t„)],  Y[y(t  )], 

/v  +4*  O  0 

Z[z,  t  ],  V[v,  t0],  V[v],  y C4> , t03 ,  and  yC^].  Here  the  arguments  such  as 
x(tQ)  are  the  statistically  contemplated  values  as  distinct  from  the  most 
probable,  or  expectation,  values  such  as  <x(tQ)>.  For  the  distributions 
such  as  Z[z,tQ],  the  quantity  z  is  time- independent,  but  the  uncertainty  in 
its  measurement  remains  time- dependent  because  it  depends  on  the  penetrator 

A 

position  relative  to  the  sensor.  Thus,  the  distribution  Z[z,tQ]  is 
time- dependent,  as  indicated  by  the  notation. 

3.  The  Distribution  Functions 


a.  Distribution  of  Measurables.  Let  us  first  enumerate  the 
distribution  functions  pertaining  to  measured  quantities.  We  take  the 
position  distribution  functions  (PDFs)  of  x(tQ)  and  y( t0 )  to  be 
Gauss i an 


X[x(t0)]  = 


1 

°x{to)V?F 


Y[y{t0)]  • 


Cy(t0)v/2ir 


exp 


exp 


•  .  [x(t0)-<x(t0)>]2  | 

'  afro)  > 

.  [y(t0)  -<y(t0)>]2  ) 

25y(t0)  I 


(1) 


(2) 


1.  The  implications  of  doing  so  are  explained  in  D.C.  Baird,  Experimentation: 
An  Introduction  to  Measurement  Theory  and  Experimental  Design.  Prentice-Hall, 
Inc.,  Englewood  Cliffs,  New  Jersey,  1962,  pp.  6-45.  (Short  Title  -  Baird) 
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We  take  a  truncated  Gaussian  for  2  [2,  t  ].  The  truncation  simulates  the 
fact  that  the  sensor  never  indicates  the  penetrator  is  at  or  below  ground 
level.  Renormalization  then  yields^ 


zG,t0] - i - — 

0  »(tn)VS 


exp 


•%s>£ 

2^lo> 


(3a) 


\  =  °2  Q(-«>/az) 


(3b) 


We  also  take  a  Gaussian  for  the  speed  distribution  function  (SDF)  of  v: 


V 


]  = 


exp 


r  ^  ^  t  2 

-  [  V-<V>] 

^(V 


(4) 


The  heading  distribution  function  (HOF)  of  p  should  reflect  the  property  that 
the  raid  enters  the  defended  area  from  a  most  favored  direction.  To  account 
for  this,  we  take  a  Gaussian  in  tan  <J): 


?[tan  *,t0]  = 


exp 


-  [tanj)  -  tan<$>]^ 
2*(tJ 

ip  0 


(5) 


1.  The  Gaussian  probability  function  Q(x)  is  tabulated  in  M.  A.  Abramowitz 
and  I. A.  Stegun,  Handbook  of  Mathematical  Functions,  National  Bureau  of 
Standards,  Washington  DC,  1965,  pp.  931-932,  966-972;  see  also  equation 
(26. 2. IS).  (Short  title  -  A&S) 
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The  assumed  HDF  is  designed  to  disallow  headings  in  the  interval  * 

Equation  (5)  also  includes  the  implicit  approximation  <tan$> = tan<$>  which  is 

A  A 

accurate  for  <1 |>>  =  0,  but  also  for  <\p>-ir. 

b.  Distribution  of  Maneuvers.  Let  us  now  enumerate  the  distribution 
functions  pertaining  to  maneuver-induced  quantities.  We  first  note  that  these 
are  likely  to  be  symmetric;  that  is,  the  penetrator  is  equally  likely  to 
change  to  a  more  westerly  heading  as  to  a  more  easterly  one,  and  it  is  equally 
likely  to  slow  down  as  to  speed  up.  These  considerations  suggest  to  account 
for  maneuvers  by  using  Gaussians  with  the  corresponding  expectation  values  set 
to  zero.  We  therefore  select  for  the  change-of-speed  distribution  function 
(CSDF ) 


vCv3  =  — i -  expffi  2 

^  2at 


(6a) 


<v>  =  0 


(6b) 


and  for  the  change-of- heading  distribution  function  (CHDF) 


+  + 


'{'[tamj;]  =  — p -  exp 


a  VZF 

<P 


-  CtamL>  -  <tantjp,]‘ 
,+2 
% 


(7a) 


<  t  anip>  =  0  (7b) 

The  CHDF  also  simulates  the  most  favored  direction  property.  In  addition,  we 

■f  ^ 

take  cfy  and  av  to  be  independent  of  one  another  and  to  be  constant  over 
the  entire  expected  path  L. 


4.  Expectation  Values  and  Standard  Deviations.  Our  enumeration  of  the 

distribution  functions  is  incomplete  without  an  estimate  of  the  expectation 

values  <x(tQ)>,  <y( tQ )> ,  and  the  maneuver- induced  standard  deviations 

a  ,  a,.  Let  us  first  show  how  to  determine  <x(t.)>,  <y(t_)>. 
v  v  oo 
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a.  Expectation  Values.  The  expectation  values  of  the 
penetrator  position  coordinates  will  be  calculated  from  the  equations 

<x(tQ)>  =  XF  +  ve  sin4»  (t0-tF)  (8) 

<y(tQ)>  =  yp  +  vg  cos<4»  (t0-tF)  (9) 

which  follow  from  figure  4.  In  these  equations,  vg  is  the  "effective 
penetrator  speed  along  L".  Even  though  the  most  probable  sensor-observed 
penetrator  velocity  <v>  lies  along  L,  its  magnitude  is  not  the  effective 
penetrator  speed  along  L;  that  is,  <v>^vg.  This  is  a  maneuver  effect  which 
has  the  consequence  that  vg  is  the  root-mean-square  of  <v>  cosj;  that  is, 


ve  3  J  cos1 2  }  exp  d(tanj) 

G\p  v'Eir  ^  > 

(10) 

.  <;>2  f  exp<-  U2<1  d  + 

'  5*  i  *  I2 


The  integral  has  a  closed-form  solution^-;  we  find 


r-—  exp(l/2o|)  +  _ _ 

- - -  erfc(l/o^  /^T) 


CTip 


<V> 


where  erfc(x)  is  the  complementary  error  function. ^ 


(ID 


1.  A&S,  equation  (7.4.11),  p.  302. 

2.  A  tabulation  of  erfc(x)  is  given  in  A&S,  pp.  297,  310-311;  see  also 

equation  (7.1.28). 
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b.  Standard  Deviations 


1.  Maneuver- Induced  Standard  Deviations.  We  next  turn  to  the 
maneuver- induced  standard  deviations  and  ay.  Even  though  we 
consider  these  parameters  as  freely  selectable,  we  can  estimate  their  values 
as  functions  of  the  penetrator  range  R  .  Figure  5  illustrates  an  extreme 
maneuver  which  takes  the  penetrator  far  away  from  the  expected  path  L.  In 
this  case,  the  travel  distance  from  Pp  to  Py  is  2a,  and  the  total  travel 
distance  is  subject  to  the  restriction 

^BF  +  +  =  Rp  (12) 

from  which 

a  =  ^(Rp  -  Rgp  ~  Rjr )  (13) 

Hence  the  angle  E,  from  figure  5  is 

S  =  arccos(RpT/2a)  (14) 

A  reasonable  estimate  of  a^is  then  given  by  the  equation 


cty  =  tan(-jS)  (15) 

Even  though  we  take  and  $v  to  be  independent,  we  can  still  derive  a 
baseline  value  for  av  through  practical  considerations  as  follows.  The 
generic  scenario  calls  for  simultaneous  weapon  detonations  which  means  that  a 
penetrator  flying  along  the  extreme  path  of  figure  5  would  have  to  move  at 
increased  speed.  We  therefore  set 

(<v>  +  oy)(tT  -  tF)  =  2a  (16) 
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Figure  5.  "Extreme  Maneuver"  -  An  Example 
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By  the  same  logic,  a  penetrator  which  flies  along  the  expected  path  L  would 
have  to  travel  at  decreased  speed.  Hence  we  write 


(<v>  -  $v)  (tT  -  tF)  =  RpT 


(17) 


The  solution  of  these  equations  is 


+ 


/2a  -  RFT\ 

\2a  +  Rft; 


A 

<V> 


(18) 


2.  Sensor-Inducea  Standard  Deviations.  We  begin  with  the 
observation  that  sensor-inGuced  standard  deviations  are  usually  expressed  in 
terms  of  ciown-range  and  cross-range  position  errors,  and  in  terms  of  radial 
ana  tangential  velocity  errors.  In  other  words,  sensor  errors  are  usually 
quoted  relative  to  a  spherical  coordinate  system  with  the  origin  at  the 
sensor.  Hence  it  is  necessary  to  transform  the  standard  deviations  from 
sensor  coordinates  to  air  base  coordinates.  The  details  are  explained  in 
appendix  A.  The  results  are: 


5x(t0)  =  ctr  sin<:[9(t0)]cos2[$(t0)]  + 

+  [aELrSp(to)]2  cos2[e(to)]cos2[$(to)]  +  (19a) 

+  taAZ  rSp(to)]2sin2[e(t0)]sin2  [4>(t0)] 

°y(t0)  =  aR  si^CQUotfsin"  [$(t0)3  + 

+  C°ELrsp(to^2cos2^^to^sin2^^to^  +  (W 

+  [aAZ  rSp(to)]2sin2[0(to)]cos2[$(to)] 
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a 2 ( tQ)  =  aR  COS2[  0(tQ)]  +  C°EL"rSP(to^2  sin2  CBit0)l  (19c) 

SJUq)  =a$R  Sin2[0  (t0)]  sin2[J  +  $(tQ)]+ 

+  av^Lcos2[0(to)]sin2[J  +  $( tQ)]  +  Cos2  [$  +  $(tQ)]  (20a) 


- 


!a2R  sin2[e ( tQ)]cos2[^  +  $(tQ)]  + 


+  ov^Lcos^[0(tn)]cos‘::[$  +  *(tft)]  + 


+  afkl  sin^[ty  +  ^t0)][/(v  cos  I|») 


2-v2 


(20b) 


Git  )  =  arcs  in) 


CXBR^o^  "  XBS^O^  +  _ 


rBPV  O'  JBSV  o' 


(21a) 


4>( tQ)  =  arcsin  j 


yBp(t0)  ' 


>£ 


xBP(to^ 


xBS^o^2  +  ^yBP^o^  "  yBS^o^2 


(21b) 


In  these  equations,  oR  is  the  range  error,  1S  the  elevation  error, 
aAZ  is  the  azimuthal  error,  avR  is  the  radial  velocity  error,  cr^  is 
the  elevation  velocity  error,  and  avki  is  the  azimuthal  velocity  error.  The 
"hat  symbol"  is  not  applied  here  because  it  is  obvious  that  these  parameters 
refer  to  the  sensor.  The  quantity  r^p(t0)  is  the  distance  between  the 
sensor  and  penetrator  at  time  tQ,  and  the  subscripts  "BS"  and  "BP"  stand  for 
"base-to-sensor"  and  "base-to-penetrator" .  (See  also  figure  A2  in  appendix 
A.)  The  angles  0(tQ),  $(tQ)  and  the  distances  rsp(t0),  r'Bp(tQ) 
are  obviously  time-dependent,  and  rBS( tQ)  is  time-dependent,  given  a 
mobile  sensor. 
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B.  EVOLUTION  OF  THE  POSITION  DISTRIBUTION  FUNCTIONS 


We  have  seen  how  to  characterize  a  sensor  reading  at  arbitrary  time  tQ. 
This  subsection  will  develop  the  PDFs  at  time  t>tQ,  given  no  additional 
sensor  readings  after  time  tQ.  The  solution  lies  in  a  series  of 
convolutions. 

Let  us  begin  with  the  "total  speed" 

v  =  v  +  v  (22) 

which  is  described  by  the  total  speed  distribution  function  (TSDF),  V[v]. 
This  function  can  be  derived  from  the  equation 

+oo 

V[v]  =  f  V[v|v]  ^[v]  dv  (23) 


where  V[v|v]  is  the  "conditional  TSDF".  It  satisfies  the  relation1 


V [v  |  v]  dv  =  V[v]  dv 


(24) 


Here,  we  observe  equation  (22)  and  the  relation  dv=dv  which  follows  from 
equation  (22)  under  the  statistical  condition  dv  =  0.  This  yields 


vCv  |  vl  =  vCv-v]  (25) 

It  then  follows  that  the  integral  of  equation  (23)  is  the  convolution  of  two 
Gaussians  which  has  the  solution 


VCv,  tQ]  = 


1 


(26a) 


1.  This  technique  for  interchanging  distribution  functions  is  common;  see  for 
example  F.S.  Hillier  and  G.J.  Lieberman,  Operations  Research,  Holden  Day,  Inc., 
San  Francisco,  1967,  p.  347. 
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where 


(to}  =  Sv  {to}  +  «v 

(26b) 

A 

=  <v> 

(26c) 

Using  the  same  techniques,  we  can  now  calculate  X[x(t)]  and  Y[y(t)].  (Our 
ignoring  maneuvers  in  the  z-direction  is  equivalent  to  setting  Z[z(t)]  *  Z[z,t0].) 
The  formal  solutions  are 


X[x(t)] 


+00  +00  +00 

J I J 

_CO  _OQ  —  00 


X[x't)|x(t0),$,  if] 


X[x(t0)]  x 


X  taniji, tQ]  $  [tan^]d[x(t0)]d(tanj)d(tanf)) 

+  oo  +oo  +00 

Y[y(t)]  =  JJJ  Y[y(t)|  y(tQ) , $  ,  J]  Y[y(t0)]  x 

-  CO  -oo  -oo 

x  ^ [ tan$, tQ]  $  [tan  J]  d[y(tQ)]  d(tan^)  d(tanj) 


(27a) 


(27b) 


Here,  we  have 


X[x(t)  lx(t0),5),J]  |dx(t)  |  =  vCv]  | dv  | 

(28a) 

Y[y(t)|y(t0),i,  $]|dy(t)j  =  V[v]|dv| 

(28b) 

These  allow  us  to  derive  the  conditional  distribution  functions  on  the  left 
side.  To  this  end,  we  first  observe  in  analogy  to  equations  (8)  and  (9)  that 


<x(t)>  =<x(tQ)>  +  vesin4>  (t-tQ) 

(29a) 

<y(t)>  =  <y(t0)>  +  vecos4>  (t-t0) 

(29b) 
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where  vg  is  the  effective  penetrator  speed  defined  by  equation  (11).  Since 
equations  (29)  contain  only  expectation  values,  in  other  words,  temporal 
constants,  they  hold  true  for  finite  time  intervals  (t-tQ).  If  we  remove 
the  bracket  operators,  we  obtain 


x(t)  =  x(tQ)  +  v  S(t-tQ)  (30a) 
y(t)  =  y(tQ)  +  v  C(t-tQ)  (30b) 
S  =  si  n((p  +  J)  (30c) 
C  =  cos  (<p  +  t)  (30d) 


However,  because  of  the  statistical  fluctuations  of  v,  ij;,  and 
(30)  hold  true  only  for  infinitesimally  small  time  intervals 
other  hand,  if  we  proceed  from  x(t),  y(t)  to  the  statistical 
dx(t),  dy(t),  the  temporal  condition  no  longer  applies.  If, 
we  observe  the  statistical  constants  listed  in  the  condition 
equations  (28),  we  obtain 


J,  equations 
(t-t0).  On  the 
uncertainties 
at  the  same  time, 
compartments  of 


dx(t)  =  (t-tQ)  S  dv  (31a) 

dy(t)  =  U-t0)  C  dv  (31b) 


We  repeat  that  these  are  not  temporal  differentials,  but  rather  these  are 
statistical  differentials. 


If  now  equations  (26),  (29),  (30),  and  (31)  are  inserted  in  equations 
(28) ,  we  obtain 

X[x(t)  |x(t  ),  J]  = 


<*-y 


exp 


[x( t )  -  x(tQ)  -  Sp(t)]‘ 

2S2avU-t0)* 


(32a) 
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Y[y(t)!y(t0),ii),ijl]  = 


(t-t0)cvlcI^T 


exp 


_[y(t)-y(t0)-cP(t)] 

T 


(32b) 


p(t)  5  <x(t)>  -  <x(tQ)> 
<S> 


<y(t)>  -  <y(t0)>  =  <v>(t_tj  (32c) 
<c>  0 


If  next  we  insert  equations  (32),  (1),  and  (2)  into  equations  (27),  we  see 
immediately  that  the  integrals  over  x(tQ)  and  y(tQ)  are  merely  simple 
convolutions  resulting  in 


X[x(t)|$,  $]=  - - -  exp 

CTS(t)  /2t 


[x(t)-<x(tQ)>  -  Sp(t)]2 
2a|(t) 


(33a) 


Y[y(t)|*,V]  = 


ac(t)  /2tt  ’ 


exp 


[y(t)-<y(t0)>-CP(t)]‘ 

2a|(t) 


a|(t)  =  a2(t0)  +  S2(t-t0)2  a2(t) 


(33b) 


(33c) 


Or(  t)  =  ^(U  +  C2(t-t)2  0j(t) 


y'  o 


For  further  treatment,  we  introduce  the  abbreviations 


A  A  ^ 

4=  t  snip  \  p  = 


We  then  note  that 


s(D.  p)  =  -A*  M 


^  (i+  t2) 


C(y,  J)  = 


I  A  + 

u 


(33d) 


(34a, b) 


J(l+uh  a+J?j 

Our  problem  now  reduces  to  solving  the  following  convolution  integrals: 


(35a, b) 


+  00  +00 


X [x ( t ) ]  =  j  J  Gx(C,t)  ?(£)  ^(S)  dydS 


_  00  .00 
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(36a) 


Y[y(t)]  =  f  j  ,  p)  ip(p)  ¥(J)dy  du 


(36b) 


+  1  [x(t)  -  <x(t  )>  ~S (y ,y ) p ( t ) ] 2 

Gx(u,  t)  =  "  A-  exP  "  - - —  (36c) 

a$(UiU)  /2tT  2a2(y,y)  ) 


G  (u,  y)  =  - - — - 

y  oc(y,u) 


Cy(t)  -  <y(t  )>  -C(u,u)p(t)]2 
exp - V— -  (36d) 

2  °2C(v,t) 


These  are  not  simple  convolutions  because  G  and  G  depend  on  y  and  y 

x  y 

in  a  "non-Gaussian"  way.  Nevertheless,  the  integrals  can  be  solved 
approximately  via  an  expansion  about  moments.  As  is  shown  in  appendix  B,  the 
solutions  to  second  order  in  and  o^  are  as  follows: 


X[x(t)]  =  Gx(4>,  o)  [1+  2Xi  (ax)1] 


(37a) 


Y[y(t)]  =  G y(4>,o)  [l-t-y^Y.(Ay)1] 


(37b) 


Ax  =  x(t)  -  <x(tQ)>  -  p(t)  sin  4> 


(37c) 


Ay  =  y(t)  -  <y(tQ)>  -  p(t)  cos<^> 


(37d) 


Xi  =  I  aS  ^  ^4xi  +  42xi^ 


Yi  ?ac  ^ 
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I 


X4  - 


aj  ( t-tQ) 4  sin2(2<i>)  cos4<J> 


4o^  ^2i f 


h  - 


cos5<^>si  n(  2<4»>) 


^7 


CTwU-02  cos4<45>( l-4si n2<ip>)  +  p2  cos®<^>  ,  ^ 

_v - £ - -  .  6ai  X, 

S  /o~“ . 1  5  4 

ac  /27T 


xi  ■ 


-3psin<^>cos  <i> 


'  °s  +  2  ‘yt-y 

.  <4  */2tt 


2  2  - 
COS  <ip> 


A  /V 


XQ  =  X4  +  X2) 


+  o}(t-t0)4  sin2(24>) 

x4  =  - 


4 /ST 


O  O 

+  a„  ( t-t  )  pcos«j»  sin  (2<^>) 
x3 - T 


X2  = 


a2  (t-tQ)2  cos  (2<$>)  +  p2  cos2<9>  o  + 


6cS  X4 


+ 

Xi 


-o  S1n<*>|  JS  *  6av  “-‘o*2  cos2<*’> 


>'°S  * 


+  2  2  +  + 
x0  “  ~  a$  X4  +  X2^ 
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(37  g ) 


(37h) 


(371) 


( 37  j ) 


(37k) 


(371) 


(37m) 


(37n) 


( 37  o ) 


( 37  p ) 


i 


+  _  (t-tQ)4  sin2(2<$>) 

Y4  4a?  JTT 


+ 

Y- 


+ 

Y, 


a2  (t*t0)2  P  sin<’^>  sin  ( 2<4») 
CTC  ^ 


p2 si n2<ip>  -  a2  (t-t  )2  cos ( 2<ip>)  + 

.5  ° - «°C  y4 


“C 


1  *  -p  coscijp 


a2  +  6  a2  (t-tn)2  sin2<J>" 


v'  ~  'o' 


+  o  2  +  + 

Y0  8  -  aC  (3aC  Y4  +  V 


( 37  v ) 


(37w) 


( 37  x ) 


(37y) 


(37z) 


A 
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V 


In  applying  equations  (37),  we  must  remember  that  now 


(t)  =  a2  (tQ)  +  sin2<^  oj(t)  ( t-tQ)2 


(t)  =  a2  (tQ)  +  cos24>  a2(t)  (t-tQ)2 


C.  MORE  ON  THE  POSITION  DISTRIBUTION  FUNCTIONS 


(37aa) 


(37ab) 


The  expansion  presented  by  equations  (37)  is  accurate  f or  pa^«as  and 
p 0^«a£.  Thus  it  has  limited  use  in  its  present  form.  However,  if  we  think 
about  the  kinematics  of  the  problem,  we  can  use  the  expansion  to  "back  into"  a 
more  useful  approximate  solution. 

Figure  6  shows  how  heading  uncertainty  enters  the  problem.  The  uncer¬ 
tainty  induces  a  contribution  p  sin(tan~4^)  to  the  x-axis  of  the  error 
basket,  and  a  contribution  p[l-cos(tan"4^)]  to  the  y-axis  of  the  error 

basket.  These  extend  a  and  a  roughly  in  the  conventional  way: 

X  y 


°x  TOT  ^2  +  P2  sin^tan'^) 


Oy  TOT  ^ac2  +  p2  [1  -  cos  (tan  "4^)]  2 


Since  we  are  treating  cases  for  which  Oy« 1,  we  note  that 
sin(tan_1a^)  = 
cos(tan-^)=  a2 


(38a) 


(38b) 


(39a) 

(39b) 


Furthermore,  whena^«l,  the  error  basket  will  continue  to  look  roughly  like  a 
Gaussian  ellipsoid  over  its  entire  evolution.  Thus,  we  are  led  to  expect 
solutions  of  the  form 


X  [x  ( t )]  % 


1 


'2*  (c2  +  M  cos24>p2  ojp 


exp 


-  (Ax)- 

x— “r 


L  2(a|  +  M  cost4>r?_cJ). 


(40a) 
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We  can  now  turn  to  the  mathematics;  i.e,  the  expansion,  to  find  out  what  M,  N 
are. 

When  we  go  to  equations  (37),  we  can  first  note  that  X4,  Y4  are  both 
always  small.  If  we  next  observe  that  for  most  interesting  problems 
av(t-t0)«p,  we  can  then  note  that  x3,  Y3  are  generally  small  compared 

to  X2,  Y2,  and  we  set 


(41a) 
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(41  i ) 


-  as  X2 


Xo  * 


2  + 

-4  h 


(41  j ) 


Yo  “ 


2  ; 
aC  Y2 


(41k) 


Yo  =  °C  Y2 


(411) 


Next,  we  expand  in  a  Taylor  series  a  Gaussian  of  the  form 


G(Ax)  = 


1 


V2rr  a2  +  K62 


exp 


L  2 (a2  +  K62)J 


(42) 


with  the  result 


G(Ax)  = 


r  -(ax)^ 

+i_ 

« 2  ’  lexJ- 

(ax)2 

-  e*p 

/2tT  cr<. 

L  J 

2  /2tt’ 

T  exP 

L  as  asJ  L 

4  J 

(43a) 


or,  more  simply, 


G(Ax)  -  G  j  1  - 


20c 


+  x 

74 

2aS 


(43b) 


where  G  is  the  "unperturbed"  Gaussian.  We  now  write  this  in  the  more 
suggestive  form 


G(Ax)  *  G 


1  *  aS  X2  +  (Ax)2  X2 


(43c) 
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With  the  help  of  equations  (37e),  (41  a) ,  (41b),  this  leads  us  to  the  identifi¬ 
cation 


/  V  ?  ^  i /N 2  4  ^  +?\ 

K6  =  p  cos  <ip>  (a^  cos  <i^>  +  a£) 


(44) 


We  therefore  reach  the  inductive  conclusion 


X[x(t)]  = 


exp 


a  x  v  tfr 


I"  -(ax)*  1  v 

t^TT  1 


Sx(<$>,0)Ax 


(45a) 


2.2  2  2^/2  4  ^  +2\ 
°x  -as  +  p  °os  <^>  (a|  cos  <ijj>  +  cfy) 


(45b) 


In  exactly  the  same  way,  we  conclude 


Y[y(t)] 


Cy  »X  TT 


exp 


20y 


+  Y1G  (4>  ,0)  Ay 


(45c) 


Qy  =  sin^<$>  (p^  cos4  <$>  +  p^) 


(45d) 


We  could  go  one  step  further  and  incorporate  the  X^,  Y^  terms.  If  we  did 
they  would  impact  the  Gaussian  through  small  additive  terms,  of  order  pa 


to  Ax  and  Ay.  Hence,  we  write  instead 


r 


X[x(t)]  ~ 


STi 7 


exp 


(46a) 


Y[y(t)]  “ 


Oy  /SlT 


exp 


hf1] 


(46b) 


These  expressions  are  accurate  for  and  for  pa|  <<CTx’  Pa|  ^y- 


D.  THE  EXPECTED  PROBABILITY  OF  KILL 


Let  PK  denote  the  conventional  kill  probability.  This  is  a  conditional 
probability,  for  it  depends  on  the  following: 

1.  The  probability  that  the  sensor  successfully  detects  and  tracks  the 
penetrator  (successful  vectoring). 

2.  The  probability  that  the  interceptor’s  on-board  sensor  detects  the 
penetrator  (successful  handover). 

3.  The  probability  that  the  weapon  launched  by  the  interceptor  success¬ 
fully  homes  on  the  penetrator,  given  penetrator  detection  (successful  homing). 

4.  The  probability  that  the  weapon  kills  the  interceptor,  given  that  it 
homes  (successful  kill). 

In  the  present  paper,  we  do  not  address  the  problem  of  assessing  these  proba¬ 
bilities,  but  rather  we  consider  them  as  inputs.  More  precisely,  we  assume  that 
all  of  the  probabilities  mentioned  above  are  consolidated  in  the  probability 

PK  =  PK  txA’  yA’  ZA’  lA  I  y<‘tA^  Z^A^ 

This  is  the  probability  that  the  weapon  which  is  launched  from  the  "attack  posi¬ 
tion"  at  the  "attack  time"  tft  will  kill  the  penetrator,  given 

that  the  penetrator,  at  the  attack  time  tA>  occupies  the  position  [x(t^), 
y(tft) ,z(tA)] .  The  present  paper  focuses  on  the  probability  that  the  penetra¬ 
tor  will  satisfy  this  condition.  This  probability  is  described  by  the  position 
distribution  functions  X[ x ( t A ) ]  and  so  forth.  With  the  aid  of  these  functions, 
we  then  calculate  the  "expected  probability  of  kill" 


+O0  +CO  +00 


<fVXA’yA,zA’V>=  f  f  f  PK^XA’yA,ZA’tAl  x(V’  y(tA^ 

*/-O0  -CO  •'-00 

z(tA)]  X  [  x ( t A ) 3  Y  C  y(tA)]  Z[z(tA)]  dx(tA)  dy(tA)  dz(tA)  (48) 


Here  we  recall  that  this  approach  accounts  for  both  sensor  errors  and  maneuver 
effects . 
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III.  THE  BASIC  ATTACK  EQUATIONS 


Any  air  intercept  problem  is  driven  by  velocity  and  geometry  considera¬ 
tions.  By  assumption,  the  geometry  of  our  problem  is  euclidean;  consequently, 
intercept  is  quantifiable  using  elementary  vector  analysis  techniques,  with 
the  resulting  solution  depending  largely  on  penetrator  and  interceptor  velo¬ 
cities.  This  section  develops  the  vector  analysis  tools  required  to  calculate 
attack  positions  and  attack  times  by  focusing  on  one  interceptor  successively 
engaging  several  penetrators.  For  this  purpose,  penetrator  maneuvers  are 
ignored;  their  effect  is  the  subject  of  section  IV.  This  section  begins  with 
the  intercept  geometry,  then  turns  to  the  impact  of  varied  interceptor  flight 
profiles,  and  finishes  by  "devectorizing"  the  attack  equations  for  computeri¬ 
zation. 


A.  INTERCEPT  GEOMETRY 


The  basic  intercept  problem  in  which  one  interceptor  successively 
attacks  several  non-maneuvering  penetrators  at  close  range  is  elementary. 
Figure  7  shows  an  overhead  view  of  the  intercept  geometry.  At  time  t£,  the 
interceptor  leaves  point  A*  at  which  it  previously  attacked  another  penetra- 
tor.  It  flies  with  velocity  Vj  along  the  transit  vector  Tj  to  attack  the 
present  penetrator  at  the  attack  point  A,  and  it  arrives  at  A  at  time  ta- 
In  general,  attack  takes  place  at  a  time  tA>tg.  For  example,  if  the  pene¬ 
trator  raid  is  tenuous,  the  interceptor  may  have  to  loiter  for  a  time  tl 
while  awaiting  the  penetrator's  arrival  at  A;  and  it  generally  will  have  to 
search  for  a  time  Tg,  depending  on  the  performance  of  the  fire  control  sys¬ 
tem.  At  the  same  time  t^  at  which  the  interceptor  leaves  point  A*,  the  pene¬ 
trator  leaves  point  Q  and  flies  with  velocity  ^  along  the  transit  vector  T  to 
point  A  where  it  is  attacked  at  time  t^.  The  points  Q,  A*,  and  A  are  repre¬ 
sented  by  vectors  r^lt^),  r and  r^,  respectively.  In  addition,  the 
C-3  network  may  carry  a  dead  time  for  data  processing  and  handover.  Let 
us  assume  that  the  vectoring  information  is  always  internally  adjusted  for  the 
dead  time  before  handover;  then  its  only  impact  is  on  the  PDFs  through  equa¬ 
tions  (37).  These  considerations  give  rise  to  the  following  equations: 


fcA  =  la  +  tL  +  TS 


(49) 
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(50) 


TI  ’  lta  -  ‘t>  VI 


T  *  (tA  -  tt)  v 


■A  ’  rQ(tt>  *  T 


rA  =  rA  +  TI 


(51) 

(52) 

(53) 


Equations  (49)  through  (53)  allow  us  to  calculate  t«  and  rA  which  are 
needed  in  equation  (48).  First  we  combine  equations  (52)  and  (53)  to  get 


■+■  -*•* 


Ti  =  rqCt*>)  +  T  - 


Inserting  equations  (50)  and  (51)  gives 

(ta  ~  t£)Vj  =  (*a  '  tjj,)  v  *  rA 

Finally,  inserting  equation  (49)  gives 

[tA  -  (V  \  +  TS}^I  =  rV^}  +  (tA  *  V  *  -  TA 
The  following  auxiliary  relations  simplify  the  notation: 


W  =  rn(tf_)  -  t£  v  -  r 


e  a  t,  *tl*ts 


Equation  (56)  now  becomes 


(54) 
(56) 

(55) 


(57) 

(58) 


(t#  -  a)  v,  ■  tA  v  +  w 


(59) 


Equation  (59)  is  the  basic  attack  equation  and  represents  the  intercept 
geometry.  Our  goal  is  to  use  it  to  find  the  attack  time  tft,  but  we  cannot 
use  it  directly  because  only  the  magnitude  and  not  the  direction  of  Vj  is 
known.  Thus,  we  first  square  equation  (59)  to  get  the  following  quadratic 
equation: 
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(60) 


(Vj  -  v2)t 2  -  2(©v2  +  v*W)tA  +  (©2v2  -  W2)  =  0 
There  are  two  solutions: 

tftlj2  =  A  4  ^A2-  b‘  (61a) 


O 

0  Vj  +  v-W 


B  = 


02  v2  -  W2 


(61b, c) 


The  corresponding  attack  vectors  follow  from  equations  (51)  and  (52): 


rAl,2  =  fy1*)  +  (tAl,2  -  *  (62> 

(If  the  attack  position  is  "past"  the  penetrator's  target,  then  the  kill  must 

not  be  counted  in  the  <N  >  total.) 

1\ 

An  understanding  of  the  nature  of  the  double  solution  arises  through  a 
simple  example.  Consider  a  penetrator  and  an  interceptor  which  are  flying  on 
a  head-on  collision  course  and  are  separated  by  a  distance  d,  as  is  depicted 
in  figure  8.  The  interceptor  has  two  choices.  The  first  is  to  fly  directly 
at  the  penetrator  and  meet  it,  and  the  corresponding  attack  time  is 


LA1 


v  +  v, 


(63a) 


The  second  is  to  fly  directly  away  from  the  penetrator  and  allow  the  penetra¬ 
tor  to  catch  up;  and  the  corresponding  attack  time  is 


t 


A2 


(63b) 
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Penetrator 


Figure  8.  Attack  Equation  Solutions  -  An  Example 


Of  course,  if  the  interceptor  sustains  a  velocity  Vj>v,  only  t ^  has  mean¬ 
ing.  In  addition,  if  the  penetrator  is  flying  directly  away  from  the  inter 
ceptor,  only  one  attack  time  exists,  namely 


(63c) 


If  in  addition  v>Vj,  then  no  attack  times  exist. 

Precisely  the  same  considerations  drive  the  scheduling  problem.  The  goal 
is  to  use  equations  (61)  to  first  calculate  attack  times  corresponding  to  the 
"direct  assault"  and  "catch-up"  mooes,  and  then  select  the  earliest  meaningful 
attack  time.  The  attack  time  is  meaningful  provided  it  satisfies  the  condi¬ 
tions 


39 


'64a,b) 


A2  -  B>0  ;  tA>  t£ 

and  a  "fuel  condition".  The  latter  follows  from  the  observation  that,  at 
attack  time  tA>  the  interceptor  has  traveled  the  total  distance 

RI(tA)s  (tA  -  tT0}  VI  (65) 

where  Tjq  is  the  take-off  time.  The  interceptor  can  return  to  the  air  base  if 
and  only  if 

1*1  (tA)  +  rA<_Rj  (66) 

B.  VARIED  INTERCEPTOR  FLIGHT  PROFILES 


The  interceptor  must  climb  to  altitude  after  take-off,  be  vectored  to  the 
first  penetrator,  and  acquire  the  penetrator  with  its  on-board  sensor.  If  the 
penetrator  flies  at  a  different  altitude,  and  the  nature  of  its  weapon  system 
or  engagement  policy  necessitates  the  interceptor  attacking  from  close  range, 
then  homing  must  occur.  The  most  interesting  case  involves  low-flying  pene- 
trators;  the  interceptor  then  must  dive  to  home  in,  attack,  and  return  to 
cruising  altitude  after  intercept.  Figure  9  depicts  this  flight  profile  by 
showing  the  interceptor  climbing  from  point  A*  and  subsequently  diving  to 
point  A.  The  climbing  (diving)  angle,  velocity,  and  time  are  yc  (yd), 
vc(vD),  and  t^(tq)»  respectively.  The  penetrator  altitude  is  h  (for  ease 
of  presentation  we  assume  that  all  penetrators  fly  at  the  same  altitude),  and 
the  interceptor's  cruising  altitude  is  h..  Figure  9  shows  explicitly  that 

-►  *  -4-  -► 

vector  Tj  is  composed  of  the  climb  vector  rc,  cruise  vector  r^,  and 
dive  vector  rQ.  We  now  investigate  how  this  flight  profile  impacts  tft. 


The  following  relations  arise  from  figures  7  and  9: 


TI 

=  rC  +  rCR  +  rD 

(67) 

*4* 

rCR 

=  (ta  -h  -  Tc  '  V  vl 

(68) 

rC 

=  Vc  Tq  cos  y^ 

(69) 
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Figure  9.  Interceptor  Flight  Profile  Against 
Low  Altitude  Penetrators  -  Side  View 


rD  = 

VD  TQ  COS  Y0 

(70) 

TC  = 

(hj  -  h)  /  vc  sin  Yq 

(71) 

TD  = 

(hj  -  h)  /  vD  sin  yq 

(72) 

Now,  r £  and 

rQ  both  have  the  same  direction  as  v^: 

-*  -► 

-► 

rC 

VI  .  -  .  VI 

=  rC  Vj  ’  r0  r0  v  j 

(73a 

Combining  equations  (69)  through  (73)  gives 
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(74a) 


-  v. 

rc  =  (hj  -  h)  cot  yc  -7- 


^  =  (hj  -  h)  cot  yd  4j- 


(74b) 


Upon  substituting  equations  (68),  (71),  (72),  and  (74)  into  equation  (67), 
we  find 


TI  =  Eta  -  t*  -  (ac  +  aD)  (hj  -  h)]v 


I 


(75a) 


,  f  =  vi  *  vc  C0S  YC  . 

c  VI  VC  s1n  YC  5 


=  VI  -  VD  cos  Yp 
D  vIvDSinYD 


(75b, c) 


By  comparing  equations  (75)  and  (50),  we  see  the  modified  attack  equation  reads 


rf  /  0  4.  A*+A 

"  (0  tLAG 


)]  = 


tA  v  +  W 


(76a) 


where 


1 


\aG  “  ^ar  +  an)  (^t  " 


C  D'  V"I 


(76b) 


The  modified  attack  time  solutions  are  those  of  equations  (61)  and  with  0  replaced  by 
.A**  A 


0  +  T 


LAG 


For  the  initiation  of  the  attack  sequence,  beginning  with  take-off  (h=0) 
and  proceeding  to  the  attack  of  the  first  penetrator,  we  set 


.0+1 


1  LAG  =  aC  hI  +  aD  ^hI  '  h) 


(77) 


1.  The  authors  admit  with  regret  that  this  symbol  suffers  from  elephantiasis. 
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Likewise,  to  finalize  the  attack  sequence,  concluding  with  the  last  attack 
followed  by  landing,  we  set 


ac  (hj  -  h)  +  aD  hj 


(78) 


where  q  denotes  the  total  number  of  penetrators  engageable  by  the  interceptor, 
that  is,  the  weapon  load. 


C.  DEVECTORING  THE  ATTACK  EQUATIONS 


The  attack  equations  must  now  be  written  in  component  form  to  allow 

2 

for  programming.  In  addition,  we  must  evaluate  terms  like  v  *W  and  W  ;  this 

can  be  done  using  the  penetrator  heading  tj>.  Let  us  begin  by  writing  v  and  W 

with  respect  to  the  basis  vectors  ev  and  e  in  figure  7: 

a  y 


v  =  v  sin^  ex  +  v  cos^  e 


rW  =  XQ ^ t£^ex  +  yQ^2^  ey 

"Nr  *  •>  * 

rA  "  XA  ex  *  *A 


Inserting  these  into  equations  (57)  gives 


(79) 

(80) 
(81) 


W  =  W  e  +  W  e 
a  a  y  y 

(82a) 

Wx  =  W“  vt*  S™*-  XA 

(82b) 

Wy  *  yq(t)l)  -  vtj,  COS  l|>  -  yj 

(82c) 

It  now  follows  that 


v  •  W 


Wy  v  sin  +  W  v  cos  ^ 
X  y 


(83) 


(84) 
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Inserting  equations  (83)  and  (84)  into  equations  (61)  allows  for  a  direct  cal¬ 
culation  of  the  attack  time  solutions.  The  component  forms  of  the  correspon¬ 
ding  attack  vectors  then  follow  from  equation  (62): 


?A1,2  =  CxQ(t,)+  (tA1>2  ‘  h)  v  sin^  K  + 


+  CyQ(t£)  +  (tAit2  -  l£)v  cos^  ey 


44 


IV.  MANEUVERS  AND  HOMING  SIMULATION 


Airbreathing  strategic  penetrators  generally  will  have  maneuver  capability 
and  will  use  this  capability  to  evade  intercept.  This  fact  makes  the  deter¬ 
ministic  vector  analysis  approach  of  the  last  section  too  simplistic  by 
itself,  yet  the  probabilistic  approach  from  section  II  adds  the  necessary 

sophistication.  That  approach  allows  for  homing  simulation.  Homing,  in  turn, 
3  3 

is  contingent  on  C  capability.  Given  the  C  capability,  one  can  in  prin¬ 
ciple  calculate  the  average  attack  times  and  attack  positions  by  averaging 
over  all  possible  intercept  geometries.  However,  several  difficulties  arise 
which  make  this  direct  approach  impractical.  First  is  the  fact  that  the  aver¬ 
aging  process  is  highly  complicated.  This  is  compounded  by  the  nature  of  the 
intercept  geometry.  As  we  saw  through  the  simple  example  depicted  in  figure 
8,  not  all  geometries  allow  for  meaningful  attack  time  solutions.  Conse¬ 
quently,  infinite  average  attack  times  are  sometimes  possible.  This  suggests 
to  truncate  the  averaging  process,  but  the  question  is  how?  In  other  words, 
we  need  a  selection  rule.  This  section  presents  an  indirect,  and  therefore 

3 

approximate,  scheme  which  couples  C  capability,  the  penetrator-interceptor 
configuration,  and  penetrator  maneuver  capability  to  (1)  the  definition  of  a 
selection  rule  for  viable  attacks,  (2)  the  calculation  of  the  corresponding 
expected  attack  times  and  positions,  and  (3)  the  calculation  of  the  correspon¬ 
ding  expected  PKs.  The  scheme  couples  the  probabilistic  approach  to  the 
intercept  geometry.  This  section  begins  with  the  best  case  --  continuous  vec¬ 
toring  information  update,  then  turns  to  the  worst  case  —  no  update,  and 
finishes  with  the  intermediate  case  --  periodic  update. 

A.  CONTINUOUS  VECTORING  UPDATE 

Continuous  update  of  vectoring  information  allows  for  homing  in  the 
purest  sense;  penetrator  maneuvers  drive  interceptor  maneuvers.  We  have 
already  described  penetrator  maneuvers  in  terms  of  a  perturbation  approach 
which  centers  on  a  straight  line  expected  path.  Consistent  with  this 
afproach,  we  now  discuss  interceptor  maneuvers  in  terms  of  an  expected  path. 
Doing  so  is  the  heart  of  our  scheme  for  incorporating  interceptor  homing. 

The  case  of  continuous  update  is  represented  by  figure  10.  It  shows 
the  penetrator  starting  at  point  < P(tj )>  which  is  the  expected  penetrator 
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Expected 
Penetrator  Path 


re 


10.  intercept  Geometry  - 


Continuous  Update 


Q  + 


position  at  time  from  the  analyst's  point  of  view.  The  point  <  P(tJl)>  for 
the  maneuvering  penetrator  now  replaces  the  point  Q  of  figure  7  for  the  non¬ 
maneuvering  penetrator.  In  parallel  to  equations  (8)  and  (9),  <P(t^)>  is 
given  by 


<x(t£)>  =  Xp  +  ve  sin<(jj>  (ta  -  tp) 

(86a) 

<y(t£)>  =  yp  +  ve  cos<ij)>  (tz  -  tF) 

(86b) 

At  the  same  time  tz,  the  interceptor  leaves  point  A*  where  it  attacked  the 
previous  penetrator.  The  maneuvering  of  the  present  penetrator  causes  the 
interceptor  to  maneuver,  too,  while  homing.  We  characterize  interceptor  mane¬ 
uvers  in  terms  of  an  expected  path  which  we  approximate  at  the  outset  to  be 
linear.  The  interceptor  then  moves  along  its  expected  path  with  speed  vIe<vI 
to  execute  attack  at  the  expected  attack  point  <A>.  Our  problem  now  reduces 
to  estimating  Vje.  In  the  perturbation  approach,  we  stated  that  a  given 
penetrator  was  equally  likely  to  deviate  from  its  expected  path  by  an  angle 
+ip  as  by  -J.  Figure  10  shows  two  such  deviations  at  point  <P(t^)>—  one 
corresponding  to  =  +arctan  and  the  other  corresponding  to  = 

-arctan  cy  We  structure  the  estimate  of  vJe  on  these  two  "root- 
mean-square"  maneuvers. 

The  attack  points  A^  and  A2  correspond  to  the  deviations  -arctan  a,j, 
and  +arctan  a^,  respectively.  They  are  calculated  by  using  the  basic  equa¬ 
tions  of  section  III.  The  attack  times  corresponding  to  A^  and  A2  follow 
from  equations  (61)  and  have  the  form 

tA(Ai,2)  *  tA  (,(,  =  4>  +  arctanaj  )  (87) 

The  corresponding  attack  vectors  then  follow  from  equation  (62).  They  have 
the  form 

rA(Ai,2)  =  <r(tZi)> +[tA(Alj2)  -  tA]  v  (88) 

where  < rUjJ >  is  the  radius  vector  to  <P(t^)>.  Finally,  the  corresponding 
interceptor  transit  vectors  follow  from  equation  (53): 
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(89) 


TU,2  =  rA(Al,2^  ‘  rA 

Our  selection  rule  then  is  to  require  that  the  attack  points  A^  and  ^ 
both  be  physically  real,  that  is,  that  they  correspond  to  meaningful  attack 
times  satisfying  the  conditions  of  equations  (64)  and  (66).  Should  the 
penetrator-i nterceptor  configuration  be  such  that  the  selection  rule  cannot  be 
satisfied,  then  the  attack  is  aborted. 


We  now  can  tie  down  <A>  based  on  figure  10.  The  angle  6  between  the  two 
interceptor  transit  vectors  is 


0  =  arccos 


(90) 


We  then  take  for  the  effective  interceptor  speed  along  its  expected  path 

vIe  ~  VI  cos  (91) 

Finally,  we  substitute  this  speed  and  the  effective  penetrator  speed  into 
equations  (61)  and  (62)  to  estimate  the  expected  attack  time  and  corresponding 
attack  position: 


tA(<A>)  =  tft  (<ij)>,  ve,  vJe) 

(92a) 

r A  ( < A> )  a  rA  [  tAC(<A  >) ,  vfi] 

(92b) 

It  now  remains  to  outline  the  corresponding  <P.,>  calculation.  Let  us 
3  N 

assume  that  the  C  network  carries  a  non-zero  dead  time  td.  The  effects 
of  dead  time  on  the  attack  calculations  depends  on  the  specific  nature  of  the 

3 

C  •  For  example,  if  the  vectoring  information  is  not  internally  adjusted  to 
compensate  for  dead  time,  then  the  resulting  attack  point  <A>  deviates  from 
the  expected  penetrator  position  <P(tA)>  at  the  attack  time  tA  by  a  factor 
of  roughly  vt^. 

The  major  impact  of  the  dead  time  is  on  the  PDFs  of  section  II.  For  one, 
the  sensor- induced  uncertainties  will  correspond  to  a  time  earlier  than  the 
attack  time: 
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a(tA)-cr  [ tA(<A>)  -  td  ] 


(93) 


In  addition,  we  have  in  the  PDFs 

p(tA)=vxTd  (94) 

This  equation  plays  the  role  of  equation  (32c)  in  the  < PK>  calculation. 
Integration  then  proceeds  according  to  the  prescription  of  equation  (48).  The 
resulting  <P^>  generally  will  be  largest  given  continuous  update;  this  is 
the  best  case. 

B.  NO  VECTORING  UPDATE 

Having  no  update  of  vectoring  information  precludes  interceptor  homing. 

In  this  case,  the  interceptor  attack  is  based  on  the  last  information 
received,  and  the  interecept  geometry  of  figure  7  represents  the  calculation 
of  the  expected  attack  position  <A>.  Let  us  assume  that  the  last  information 
was  received  at  time  t^  t^.  Then  the  point  <P(t^)>  in  figure  7  is  defined 

by 

<xttj,)>  =  Xp  +  [(t0  -  tp)ve  +  (ta  -  t0)<v>]  sin<i>  (95a) 

<ylt4)>  *  yF  +  [( t0  -  tp)ve  +  (t4  -  t0)<v>]  cos4>  (95b) 

The  attack  times  and  attack  positions  now  follow  directly  from  equations  (61) 
and  (62). 

The  < PK>  corresponding  to  no  vectoring  update  will  generally  be  smallest 
because  the  errors  in  the  sensor-predicted  penetrator  position  at  time  tA 
will  generally  be  largest.  This  is  the  worst  case. 

C.  PERIODIC  VECTORING  UPDATE 

Let  us  now  suppose  vectoring  information  is  periodically  received  with  the 
period  t  •  Then  the  average  separation  of  the  attack  time  and  the  last 
update  is  t  /2.  An  approximate  solution  of  the  corresponding  homing 
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problem  involves  a  combination  of  the  techniques  explained  in  subsections  A 
and  S  and  is  represented  by  figure  11.  We  first  calculate  the  attack  time 
t A  and  attack  position  A1  corresponding  to  continuous  update  from  equations 
(92).  We  next  subtract  tv/2  from  tA  to  find  the  "intermediate  time" 

tIN  E  '  V2 

At  times  between  t^  and  t^,  sensor  information  is  unavailable.  We  next 
find  the  corresponding  intermediate  attack  point  which  is  the 

expected  interceptor  position  at  the  time  tj^  from  the  analyst's  point  of 
view.  This  point,  shewn  in  figure  11,  is  given  by 


'lN>  ‘  (tIN  ’  ^  vIe  <r^> 


Here,  <rIN>  is  the  radius  vector  to  <pi(tIN)>»  and  <r^>  is  the  radius 
vector  to  <A’>.  While  the  interceptor  flies  to<Pj(tjN)>,  the  penetrator 
flies  to  the  intermediate  position  <P(tjN)>  represented  by 


<x(trN)>  =<x(tj>  +  v  sin<£>  (trw  -  t£) 


IN  "  ^ 

(98a) 

IN  '  V 

(98b) 

<y(tIN)>  =  <y(t&)>  +  ve  cos<^>  (tIN  -  tA)  (98b) 

The  points  <pj(tjfg)>  and  now  serve  as  the  initial  interceptor 

point  and  initial  penetrator  point.  Because  tjN  marks  the  termination  of 
vectoring  information,  we  now  revert  to  the  no  update  formalism.  The  attack 
time  tft  and  attack  position  <A>  estimates  then  follow  from  equations  (61) 
and  (62). 

If  is  small,  then  it  follows  that  the  major  impact  is  on  the  PDFs  in 
parallel  to  equations  (93)  and  (94),  except  now  is  replaced  with 
+  Ty/2.  The  resultant  <PK>  is  generally  smaller  than  that  of  con¬ 
tinuous  update,  but  is  generally  larger  than  that  of  no  update.  Periodic 
update  is  the  intermediate  case. 
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gure  11.  Intercept  Geometry  -  Periodic  Update 


V.  GUIDED  AAM  STANDOFF  ATTACK  —  AN  EXAMPLE 


Our  air  defense  formalism  can  be  applied  to  virtually  any  air  defense 
system  predicated  on  airborne  platforms.  To  illustrate  one  such  application, 
let  us  consider  standoff  attack  against  the  strategic  penetrators.  Specifi¬ 
cally,  let  us  consider  a  hypothetical  futuristic  weapon  system  based  on  long- 
range  air-to-air  missiles.  The  weapon  platform  is  given  continuous  vectoring 
information  update.  This  section  outlines  the  application  of  our  formalism  to 
model  this  air  defense  system. 

The  missile  flight  makes  the  representation  of  standoff  missile  attack 

slightly  more  complicated  than  that  of  close-range  attack.  In  figure  12, 

which  shows  the  standoff  attack  intercept  geometry,  the  attack  point  A  is 

detached  from  the  penetrator's  path.  The  missile  flies  from  A  along  vector  m 

to  kill  the  penetrator  at  the  kill  time  t^  and  the  kill  position  K  which  is 

represented  by  the  vector  r^.  The  magnitude  of  m  is  taken  as  a  parameter 

which  depends  on  the  operational  missile  firing  range  R.  (generally  less 

M 

than  the  missile's  total  range  capability): 

m  =  Rm  (99) 

Likewise,  in  the  absence  of  penetrator  maneuvers,  the  missile  flight  time  is  a 
parameter  depending  on  R^  and  the  missile  speed  vM, 


rM  =  Ru/v 
M  MM 


(100) 


Apart  from  these  missi le-inauced  modifications,  the  intercept  geometry  is  the 
same  as  that  for  close-range  attack  (see  figure  7). 


Given  that  the  penetrator  maneuvers,  the  present  problem  decomposes  into 
two  separate  homing  problems.  The  first  is  the  homing  of  the  interceptor 
toward  the  expected  attack  position  <A>.  The  second  is  the  homing  of  the 
missile  f rom  < A  > toward  the  expected  kill  position  <K>.  We  begin  with  the 
calculation  of  <A>  in  parallel  to  sections  III  and  IV. 


In  parallel  to  equations  (49)  -  (53),  the  intercept  geometry  gives  rise 
to  the  following  relations: 
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t*  ■  ‘a  *  tm  *  tl  *  Ts 


(101) 


‘a  ‘  *K  -  XN  <102> 

T*  ■  (t,  -  tt)  V*  (103) 

T  =  (fcK  -  t&)  v"  (104) 

"K  =  ^Qlt4)  +  T  (105) 

rK  =  r"A*  +  Tj  +  lit  (106) 

To  calculate  the  attack  time  tft,  we  combine  these  equations  and  find 

(tA  -  0)  Vi  +  in  =  tA  v  +  WM  (107a) 

WM  5  rQ(t4)  -  (te  -  tm)  v  -  rA*  (107b) 

0  5  h  +  TL  +  Ts  (107c) 

Squaring  equation  (107a)  gives 


O  r 


2(9v j  +  V-  WM  -  VjRM)tA  +  (rjj-j  +  q2m\  +  2evzRM  -  wj)  ■  0  (108) 


The  solutions  are 


tAl ,  2 


=  AM± 


+  v  •  WM  -  VjRm 


(109a) 


v 


(109b) 


These  solutions  do  not  reflect  penetrator  maneuvers.  Rather,  they  are 
tools  to  be  used  in  the  homing  calculations  set  forth  in  section  IV.  The 
expectea  attack  time  tA(< A>)  and  attack  position  rft(<A>)  follow  from  the 
recipe  represented  by  equations  (92).  The  expected  attack  position  <A>  now 
forms  the  basis  for  calculating  the  expected  kill  position  <K>.  The  problem 
geometry  corresponding  to  missile  flight  is  identical  to  the  close  range 
intercept  geometry  of  figure  7.  Thus,  we  merely  replace  A  by<A>,  Vj  by 
vM,  and  t£  by  tA(<A>),  and  then  use  the  formalism  of  sections  III  and  IV 
directly  to  calculate  <K  > and  tK ( < K> ) .  The  expected  attack  times  and  posi¬ 

tions  impact  <PK>  through  equation  (48)  and  interceptor  scheduling,  while 
the  expected  kill  times  and  positions  impact  <NK>. 
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VI.  SUMMARY 


The  methodology  presented  in  this  paper  is  designed  to  model  a  complicated 
strategic  defense  configuration.  It  relies  on  a  wide  variety  of  inputs  to 
characterize  geography,  force  structure,  and  system  parameters.  For  modelling 
the  close-range  (one-on-one)  intercept  problem  the  required  inputs  are: 

Geography’'’  -  fence  configuration,  coordinates  of:  air  bases, 
home/staging  bases,  US  targets,  recovery  bases,  surveillance  sensors. 

Force  Structure  -  penetrator  scheduling  against  US  targets,  number 
of  penetrators  per  home/staging  base,  number  of  interceptors  per 
interceptor  base. 

System  Parameters  - 

penetrator:  Rp,  v,  h, 

interceptor:  Rj,  Vj,  hj,  Yc,  Yp,  Vp,  Vp,  q,  Pk  envelope 
surveillance  sensor:  aR,  aEL,  CTAZ,  avR,  avEL,  avAZ,  Pp  envelope 


Some  of  the  inputs  may  require  separate  models;  examples  are  PK  and  x^. 

The  calculations  of  expected  attack  times  and  positions  and  expected  kill 

2 

probabilities  now  proceed  along  the  following  sequence  : 

Step  1  -  Calculate  vQ  from  equation  (11). 


i.  Depending  on  user  needs,  all  geography  inputs  can  be  either  abstract  or 
derived  directly  from  force  structure  and  weapon  laydown  considerations. 

For  concreteness  we  have  set  forth  the  sequence  which  corresponds  to  the 
continuous  vectoring  update  case 
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Step  2  -  For  each  interceptor  -  penetrator  configuration,  calculate 
the  earliest  tA(A)  from  equations  (87)  and  (61).  Abort  attacks 
not  satisfying  the  conditions  set  forth  in  equations  (64)  and  (66). 

Step  3  -  For  interceptor  -  penetrator  configurations  surviving  Step 
2,  calculate  rA(A)  from  equation  (88). 

Step  4  -  Calculate  Tq  ^  from  equation  (89). 

Step  5  -  Calculate  v^  from  equation  (91). 

Step  6  -  Calculate  the  earliest  t^(<A>)  from  equations  (92)  and 
(61).  Abort  attacks  not  satisfying  the  conditions  set  forth  in 
equations  (64)  and  (66). 

Step  7  -  Calculate  rA(<A^  from  equation  (62). 

Step  8  -  Calculate  ox,  3 i  3 v  Sy,  and  from  equations  (19), 

(20),  and  (21). 

Step  9  -  Calculate  < P« ( x/y»  yA,  zA)>  from  equation  (48). 

In  summary,  the  methodology  combines  surveillance  sensor,  interceptor,  and 

3 

C  capabilities  into  an  overall  air  defense  simulation  model  to  serve  as  a 
basis  for  assessing  performance  of  weapon  systems  relying  on  airborne  plat¬ 
forms.  Because  it  can  accept  any  geographic  configuration,  be  it  abstract  or 
real,  it  allows  for  modelling  a  variety  of  scenarios.  Furthermore,  it  incor¬ 
porates  maneuver  capability  of  airbreathing  penetrators  into  a  straightfor¬ 
ward  probabilistic  approach  to  generate  closed  form  solutions  which  are 
amenable  to  direct  scheduling.  In  this  way,  the  methodology  facilitates 
efficient,  time-saving  programming.  Finally,  and  perhaps  most  importantly,  it 
provides  the  tools  for  a  direct  sensitivity  analysis  of  overall  weapon  system 
performance  to  tracking  errors  and  other  determinants. 
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APPENDIX  A 


TRANSFORMATION  OF  STANDARD  DEVIATIONS 


This  appendix  details  the  calculations  which  transform  sensor-induced 
aown-range  and  cross-range  standard  deviations  into  the  cartesian  standard 
deviations  listed  in  equations  (19),  (20),  and  (21). 

A.  CALCULATION  OF  ax{tQ),  Sy(tQ),  az(tQ) 

Figure  Al  shows  the  cartesian  coordinate  system  (x^,  y^,  z^)  centered 
on  the  sensor  S.  The  penetrator  is  at  point  P(t  )  with  the  coordinates 
f>sp(t0),  ysp(t0),  zsp ( t0 ) » 1  -  The  following  relations  derive  from 
the  figure: 

xsp(t0)  =  rsp(t0)  s1n[e(t0)]  cos[*(t0)] 

(Ala) 

ySp(t0)  =  rSP(t0)  sinC<t>(t0)3 

(Alb) 

2SP(t0)  =  rSP(V  cos^^o^ 

(Ale) 

One  can  show  that  the  standard  deviation  associated  with 

is  itself  a  function  of  many  (n)  variables  can  be  calculated 
a  summation  process  as  follows*: 

a  quantity  m  which 

approximately  via 

5«  ■  s»? 

(A2) 

where  the  coefficients  of  the  standard  deviations  aYi  are  evaluated  by  sub 

A  1 

stituting  measured  values  (hence,  the  presence  of  the  hat-symbol).  Direct 
application  of  equation  (A2)  gives 


1.  Baird,  pp.  61-64. 
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1 

l(*Ks  f  , 

(**s  ) 

2 

? 

/9XS 

f  2 

axS  =  » 

\3rSP/  °R  + 

\  39  / 

°EL  + 

/  aAZ 

(A3) 


and  likewise  for  a ^  and  a^.  We  drop  the  hat-symbol  from  a^,  a^L,  a ^ 
because  these  obviously  refer  to  the  sensor. 

It  remains  to  convert  to  "air  base"  coordinates.  Figure  A2  implies  that 

rSP^V  =  rBP^V  ‘  rBS^V  ^ 

where  the  subscripts  "SP",  “BP",  and  "BS"  stand  for  "sensor-to-penetrator" , 
"base-to-penetrator",  and  "base-to-sensor".  Note  that  really  is  time- 
aependent  if  the  sensor  is  mobile.  If  we  ignore  errors  in  rg<-,  then  it 
follows  for  equation  (A4)  and  figure  A2  that 

VV  =  5xS(to};  fyV  =  5yS(to);  VV  =  5zS(to}  (A5a,b,c) 

Thus,  a  direct  substitution  of  equations  (Al)  into  equation  (A3)  results  in 
equations  (i9),  and  equations  (21)  follow  from  equations  (Al)  and  (A4). 


B.  CALCULATION  OF  3 v(tQ) 


We  are  considering  penetrators  constrained  to  a  horizontal  plane;  the  car¬ 
tesian  velocity  components  are 


v  =  v  sin  tji ;  v  =  v  cos  ip 

X  j 


(A6a,  b) 


The  corresponding  velocity  components  vr(t0),  vQ ( tQ ) ,  ( tQ )  in  a 

spherical  coordinate  system  centered  on  the  sensor  are  related  to  vx  and 
Vy  through  an  orthonormal  transformation  matrix  M(tQ): 


- 

-  - 

r  -| 

-  - 

V 

V 

V 

v„ 

X 

r 

r 

T 

X 

vy 

=  M 

v9 

9 

* 

»  mt 

vy 

0 

* 

0 

-  - 

-  - 

•  J 

(A7a,b) 
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Figure  A2.  Sensor  Coordinates  Vs  Air  Base  Coordinates 
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M(t0) 


sin[e(t0)]cos[4>(t0)] 

sin[0(to)]sin[<ji(to)] 

cos[e(t0)3 


cos[9(t0)]cos[(j)(t0)]  -sin[0(to)j 
cos[e(t0)]sin[<t)(t0)]  cos[$(t0)j 


(A7c) 


-sin[e(t0)  ]  0  J 


Now,  the  penetrator  speed  is  simply 


v 


An  application  of  equation  (A2)  gives 


qv(to)gVl~4'°))  avR  + 


(A8) 


(A9) 


Substituting  equations  (A7b)  and  (A6)  into  equation  (A9)  gives  equation  (20a). 
C.  CALCULATION  OF  a  (tQ) 

Equations  (A6)  imply  that 


tan  p  = 

y 


from  which 


3  tan  4>  3  tan  ip  3  vx  3  tan  p  dvy 

dvr  =  ~  +  »r 

Using  equation  (A7a)  we  find 


/V  A  /N 

3  tan  .4;  sin[9(t0)]  cosfty  +  <j>(t  )  ] 

-  =  - - - g- - 

3Vr  V  cos  <P 


(A10) 


(All) 


(A12a) 
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In  the  same  way  we  also  find 


3t  an  ^ 


3V0 

Stamp 


cos[0(to)]  cos[<P  +Mt0)] 

“a  O 

V  COS  ip 

-sin[ip  +  $(t0)  ] 

- 7 - p - 

V  COS  ip 


(A12b) 


(Al2c) 


Direct  application  of  equation  (A2)  followed  by  substitution  of  equations 
(A12)  yielas  equation  (20b). 
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APPENDIX  B 


CONVOLUTION  INTEGRALS  AND  MOMENT  EXPANSIONS 

This  appendix  explains  the  approximation  procedure  which  converts  equa¬ 
tions  (36)  into  equations  (37).  To  this  end,  we  observe  that  equations  (37) 

/v  + 

describe  an  averaging  process  overy  and  y.  Hence  we  may  write 


X[x(t)] 

=  <Gx(y,p)> 

(Bla) 

Y[y(t)] 

II 

A 

CD 

TT  > 

V 

(Bib) 

We  first  develop  the  functions  to  be  averaged,  G  and  G  ,  as  Taylor 

~  ^  +  *  y 

series  in  y  and  y  in  the  neighborhood  of  <y>  and  <y>.  If  we  then  carry 
out  the  averaging  process  and  if  we  observe  the  statistical  independence  of 
and  y,  we  obtain 


A  X  ^  Ct  ^  y  ^  X 

X[x(t)]  =  Gx (<u> ,  <y>)  +  <Ay  >  * — (<y>,<y>)  + 

ar  2C 

„+<’UX,/'  +  ,  ~  +3  X/^  +  . 

+  <A  y  >  — —  (<y>  ,  <y>)  +  <Ay><Ajj>  V^+  ,<?>)  + 

3y  3y  3y 

+  \  <(Au)Z>  — (<y>,<5>)  + 

3y 


1 

1 


<(a  tv 


+r 

3y 


(<y >,  <y>)  + 


(B2a) 


Ay  =  u  -  <y> 


+  +  + 

Ay  =  y  -<y > 


(B2b,c) 


Here,  we  note,  for  example,  that  <  (y-<y  >)  >  =  0,  that<y>  =  0,  and  that 

A  A  2  aO 

<(y-<u>)  >  =  cr£.  If  we  neglect  terms  higher  than  second  order,  we  obtain 
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-c> 


(B3a) 


2g  s2g 

X[x(t)]  =  Gx(<4>,  0)  +  \  (<4>.  o)  +  2  a2  — ^2  (<y>,  0) 


In  exactly  the  same  way,  we  also  find 


1  ~2  3  Gy  ,  -  ^  1  +2  3  Gv 


Y[y(t)]  -  G  0)  +  (<u>,  0)  +  2  a\p  +2  ( B3b ) 


It  remains  to  calculate  the  derivatives.  To  this  end,  we  write  equation 
36c )  in  the  shorthand  form 


JTi r  Gx  =  a  e' 


A  x2(t) 

^as  lt> 


Ax(t) =  x(t)  -  <x(t0)>  -  Sp(t)  (B4d) 

^  + 

We  denote  derivatives  with  respect  to  both  v  and  m  by  primes.  Differentiating 
equation  (B4a)  then  yields 


/2tF  g"  =  la"  -  2a ' b '  -  ab"  +  a  b‘2)  e‘b 


From  equations  ( 84b }  and(B4c)  we  get 


<2  » 

-  2a4  -aa 

a - 7 - 


(66b) 


o 
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b  = 


-  .  2 
o A x  A x  -  a  Ax 


b"  =  ~  j_qgAX^L  +  g2AX  2  +^2  Ax  A x“  -gg"AX2  +  3g2  AX2 

a 


Equation  (B4d)  yielas 


(B6c) 


(B6b) 


A  X 

=  -  s  p 

(B7a) 

it 

Ax 

=  -  s"  p 

(B7b) 

Finally,  from  equation  (33c),  we  obtain 

1 

a2  (t-tj2  S  S' 

a  = 
S 

V  0 

as 

(B8a) 

a 

a  - 

t  i2r  2  .  2  c ' 2  2/+  +  n2  c2  c  1 2-. 

^  +a  s  S  “Cv(  t-tQ)  S  S  J 

(B8b) 

S 

3 

aS 

If  now  equations  (B6),  (B7),  and  (B8)  are  substituted  into  equation  (B5),  it 
follows  after  some  calculation  that 


1  2  2  2  2  1 2  11 
(posS  )  +  °v  (t-tQ)  °s  (S  +  SS  ) 

X2  =  07  ^ 

paj  S"  -  6pav2(t-t0)2  S  s'2 
xi  =  - - 

as  /2it 

Xo  =  -  ° $  (3as  X4  +  X2) 


6a^  X4 


(B9d) 


(B9e) 


(B9f) 


For  the  evaluation  of  the  functions  X^,  it  is  necessary  to  calculate  the 
derivatives  of  S.  From  equations  (35a)  and  (35b),  we  obtain 

S(<u>,  0)  =  sin  <ip>  (BlOa) 

■—  (>y>  ,  0)  =  cos3<$>  ;  ^§-  (<U>,  0)  =  cos  <$>  (B10b,c) 

3u  3  u 


2 

(<  u>,  0)  c  -3  sin <$>  cos*  <5,> 

3S2 

2 

(<u>,  o)  =  -sin <<ij> 


(BlOd) 


(BlOe) 


Now,  returning  to  equations  (B9),  we  note  that  it  will  be  necessary  to 
distinguish  between  differentiations  with  respect  to  u  and  u  •  We  achieve 
this  by  means  of  the  hat  and  the  cross  symbols.  We  have,  for  example. 


~ll 


Xi  (Ax)1 


(Bll) 


Here  the  functions  Xi  are  obtained  by  substitution  of  equations  (BlOa), 

(BlOb)^  and  (BlOd)  into  equations  (B9).  In  a  similar  way,  we  obtain  the  func¬ 
tions  X..  The  results  are  represented  by  equations  (37g)  through  (37p). 
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An  identical  set  of  calculations  yields  the  approximat 

A  +  *  + 

functions  corresponding  to  the  and  X.  are  now  Yi  and  Yi 
represented  by  equations  (37q)  through  (37z). 


If  now  we  combine  equations  (B3a)  and(B5),  we  obtain 
X[x(t)]  =  Gx  (<y>,  0)  + 


+ 


1 

2 


exp 


+2 


(Ax)1' 


Here  we  introduce 


xi  *  o  aS  ^  Xi  +  xi ) 


Equation  (B3a)  may  then  be  written  in  the  form 


X[x(t)]  =  G  (<y>,  0) 


4 

i  +  Xi  (AX)1 


It  remains  to  show  that  equations  (37)  are  normalized, 
introduce  the  normalization  constants 


/ 


Nx  ■  |  X [x ( t ) ]  dx  (t) 


+00 


Ny  =  I  Y[y(t)]  dy  (t) 


ion  of  Y[y(t) ].  The 
.  They  are 


(B12) 


( B13 ) 


(B14) 


To  this  end,  we 


(Bl5a) 


(B15b) 
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We  then  have 


Gx  (<u>,  0)  dx  +  ir 


(<u>,  0)  dx 


(B16) 


and  a  corresponding  equation  for  Ny.  If  we  next  contemplate  substitution  of 
equations  (B9)  into  equation  (B16),  we  get  five  types  of  integrals  in  dx  or, 
what  amounts  to  the  same,  in  d(Ax).  These  are 


+oo 


9 

i 


(Ax)3 


4 

(Ax)  exp 


_  00 


/  AX2\. 


d(Ax)  = 


(BI7a) 


(B17b) 


(B17c) 


(B17d) 


(B17e) 


Hence  it  follows  that 


Koo 

J 


Gx  (<u>>  0)  dx 


(B18) 
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I  If  equations  (37g)  through  (37k)  are  observed  here,  then  it  follows  that  this 

i  +n 

integral  vanishes.  The  same  results  arise  upon  integrating  over  Gx* 

Hence  we  have 

Nx  =  1  (B19a) 

II 

The  same  operations  on  Gy  give 

Ny  =  1  (B19b) 


APPENDIX  C 


THE  APVAST  COMPUTER  PROGRAM 

A.  CAPABILITIES.  The  APVAST  program  is  designed  to  calculate  expected 
number  of  penetrators  killed.  It  begins  by  estimating  penetrator  maneuver 
tendency  based  on  fuel  constraints.  Then,  for  an  individual  attack,  it  incor¬ 
porates  the  penetrator  maneuver  tendency  into  the  calculation  of  the  earliest 
expected  attack  time  and  the  corresponding  attack  position  of  the  inter¬ 
ceptor.  Once  this  is  done,  the  program  turns  to  finding  the  associated 
expected  probability  of  kill.  The  calculations  incorporate  the  capabilities 

of  both  the  interceptor  and  the  interceptor  weapon.  Independent  of  the  nature 
3 

of  the  C  ,  the  program  always  allows  for  interceptor  homing  in  calculating 
attack  times.  Because  this  approach  is  strictly  appropriate  only  when  the 
provides  continuous  vectoring  information  update,  there  will  generally  be 
an  error  in  the  attack  time  estimates,  given  only  periodic  or  no  vectoring 
information  update.  Nevertheless,  for  most  cases  of  interest,  the  error  will 
be  small.  Aside  from  the  attack  time  calculation,  the  program  estimates 

3 

expected  kill  probability  in  a  manner  consistent  with  the  nature  of  the  C  • 

It  incorporates  a  simplistic  "cookie-cutter'1  model  for  both  the  interceptor's 
lethal  envelope  and  the  sensor  coverage.  For  example,  the  conventional  kill 
probability  is  taken  to  be  zero  outside  the  weapon  range  and  constant  inside 
the  weapon  range,  independent  of  azimuthal  angle.  At  his  discretion,  the  user 
can  fold  in  more  complicated  envelupes. 

Penetrators  are  assumed  to  fly  at  constant  altitude  from  the  point  at 
which  they  penetrate  the  warning  network  fence  until  they  reach  their  target 
location.  The  amount  of  maneuvering  that  the  penetrator  does  is  determined  by 
the  ratio  of  the  total  amount  of  fuel  available  for  this  segment  of  the  flight 
to  the  amount  needed  to  fly  the  segment  in  a  straight  line.  This  ratio  should 
be  greater  than  or  equal  to  one,  or  the  penetrator  would  run  out  of  fuel 
before  it  reached  its  target.  Since  the  degree  of  maneuvering  determines  the 
effective  velocity  of  the  penetrator  between  the  fence  penetration  point  and 
its  target,  .he  program  calculates  the  fence  penetration  time  based  on  the 
time  the  penetrator  is  to  be  at  the  target.  The  program  also  has  the 
capability  to  simulate  simultaneous  target  attacks  by  all  the  penetrators. 
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The  program  determines  its  own  start  time  for  the  scenario.  The  start  time  is 
the  earliest  penetrator  fence  crossing  time  plus  a  scramble  time  for  the 
interceptors.  No  penetrators  are  attacked  after  they  reach  their  target. 


Although  the  coordinate  system  for  the  model  is  cartesian,  the  tracking 
sensor  coverage  is  simulated  in  a  spherical  coordinate  system  to  properly 
account  for  horizon  masking  of  low  penetrators. 

Three  types  of  interceptor  bases  may  be  specified  in  the  model:  (1)  bases 
with  refueling  and  weapon  loading  facilities,  (2)  bases  with  refueling  facili¬ 
ties,  and  (3)  bases  with  runways  only.  The  model  will  always  insure  that  the 
interceptor  is  able  to  land  at  a  base  upon  completing  its  mission.  The  model 
will  try  to  find  a  base  for  the  interceptor  to  land  in  the  order  given  above. 
If  an  interceptor  lands  at  a  type  1  base,  it  is  refueled,  reweaponed,  and  is 
again  ready  to  engage  penetrators.  If  an  interceptor  with  no  remaining  wea¬ 
pons  lands  at  a  type  2  base  and  there  is  a  type  1  base  within  its  operating 
range,  the  interceptor  is  flown  to  that  base,  refueled,  reweaponed,  and  is 
again  ready  to  engage  penetrators.  Otherwise,  it  is  dropped  from  the  prob¬ 
lem.  If  the  interceptor  has  some  weapons  left,  it  is  refueled  and  is  again 
ready  to  engage  penetrators. 

If  the  interceptor  is  forced  to  land  at  a  type  3  base,  it  is  dropped  from 
the  problem  because  of  lack  of  fuel.  At  the  start  of  the  problem  the  inter¬ 
ceptors  are  assumed  to  be  fully  loaded  with  fuel  and  weapons,  and  can  there¬ 
fore  use  any  type  of  base  initially.  The  model  only  accommodates  interceptors 
initially  located  at  one  base  in  a  single  problem.  Problems  with  multiple 
basing  are  easily  accommodated  by  running  multiple  problems  that  build  on  the 
results  of  previous  cases.  This  is  done  by  not  zeroing  the  kill  probabilities 
for  the  penetrators  between  problems  and  is  user  controllable.  The  first 
interceptor  is  able  to  takeoff  at  the  problem  start  time, and  successive  inter¬ 
ceptors  at  specified  intervals  thereafter.  The  model  assumes  there  will  be  no 
landing  and  takeoff  conflicts  between  the  interceptors  later  in  the  problem. 
The  interceptors  will  not  actually  leave  the  base  for  their  first  engagement 
until  the  target  first  enters  a  tracking  sensor  coverage  so  that  the  inter¬ 
ceptor  can  be  vectored  to  the  engagement  point.  The  interceptor  is  always 
vectored  to  the  target  that  can  be  engaged  earliest,  but  it  cannot  attack  the 
same  target  successively.  Targets  are  considered  killed  when  their  cumulative 


74 


kill  probability  exceeds  a  user-determined  confidence  level.  The  interceptor 
can  attack  from  it  cruise  altitude  or  can  be  forced  to  dive  to  the  penetrator 
altitude  for  the  attack  and  return  to  the  cruise  altitude  between  attacks. 

B.  INPUT.  Input  to  the  APVAST  program  is  based  on  cards  containing  an 
identifier  and,  generally,  data  associated  with  the  identifier.  The  identi¬ 
fier  may  be  the  name  of  a  variable  in  the  model.  The  first  ten  spaces  on  the 
card  are  used  for  the  identifier,  left  justified,  followed  by  two  ten  column 
fields  for  floating  point  data.  The  remainder  of  the  card  may  be  used  for 
comments.  In  many  instances,  this  format  is  not  adequate  to  allow  the  speci¬ 
fication  of  all  data  associated  with  the  identifier.  In  these  cases,  the 
identifier  card  is  immediately  followed  by  additional  specially  formatted  data 
cards.  The  order  of  the  identifier  cards  within  a  data  deck  is  not  important, 
although  cards  associated  with  an  identifier  card  must  be  in  proper  order.  A 
specific  case  is  set  up  by  defining  all  the  necessary  information  for  the 
model  with  the  identifier  cards  and  their  associated  cards.  The  case  is  ter¬ 
minated  by  the  identifier  "ENDCASE."  Multiple  cases  may  be  run  by  the  model. 
Once  data  are  defined  by  an  identifier  card,  they  may  be  changed  by  specifying 
new  data  with  another  identifier  card  with  the  same  identifier  name.  Data 
which  does  not  change  from  case  to  case  need  only  be  input  once.  Additional 
cases  may  be  defined  by  changing  parameters  of  the  previous  data  set,  thus 
simplifying  the  generation  of  multiple  cases.  Program  execution  is  immedi¬ 
ately  terminated  by  the  identifier  "ENDJOB."  The  program  prints  out  all  the 
input  data  to  provide  a  permanent  record  of  the  parameters  associated  with 
each  run.  Data  appearing  on  cards  associated  with  the  identifier  card  are 
labeled  in  the  output  with  variable  names  or  descriptions.  Parameter  names 
and  descriptions  are  given  in  this  section  in  addition  to  the  identifier  names 
to  help  describe  the  input  to  the  model.  The  model  checks  for  invalid  identi¬ 
fier  names  and  will  issue  a  diagnostic  message  and  terminate  execution  when 
one  is  found. 

Organizations  desiring  copies  of  the  computer  code  should  complete  the 
attached  Statement  of  Terms  and  Conditions,  Release  of  Air  Force-Owned  or 
Developed  Computer  Software,  and  send  the  request  to  the  Directorate  of 
Aerospace  Studies,  Kirtland  AFB,  New  Mexico  87117. 
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STATEMENT  OF  TERMS  AND  CONDITIONS 


RELEASE  OF  AIR  FORCE-OWNED  OR  DEVELOPED  COMPUTER  SOFTWARE 

DATE 


1.  In  accordance  with  the  provisions  of  AFR  300-6,  release  of  the  following 
U.S.  Air  Force  software  package  (computer  programs,  systems  descriptions, 
and/or  documentation)  is  requested: 


2.  The  requested  software  package  will  be  used  for  the  following  purposes: 


Such  use  is  projected  to  accrue  benefit  to  the  Government  as  follows: 


3.  I/we  will  be  responsible  for  assuring  that  the  software  that  we  received 
will  not  be  used  for  any  purpose  other  than  shown  in  paragraph  2  above;  also, 
it  will  not  be  released  to  anyone  without  the  prior  approval  of  the  Air  Force. 
Further,  the  release  of  the  requested  software  package  will  not  result  in 
competition  with  other  software  packages  offered  by  commercial  firms. 


4.  I/we  guarantee  that  the  provided  software  package,  and/or  any  modified 
version  thereof,  will  not  be  published  for  profit  or  in  any  manner  offered 
for  sale  to  the  government;  it  will  not  be  sold  or  given  to  any  other 
activity  or  firm,  without  the  prior  written  approval  of  the  Air  Force. 

If  this  software  is  modified  or  enhanced  using  government  funds,  the  Government 
owns  the  results,  whether  the  software  is  the  basis  of,  or  incidental  to  a 
contract.  The  Government  shall  not  pay  a  second  time  for  this  software  or 
the  enhanced/modified  version  thereof.  The  package  may  be  used  in  contract 
with  the  Government  but  no  charge  may  be  made  for  its  use. 
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5.  The  U.S.  Air  Force  is  neither  liable  nor  responsible  for  maintenance, 
updating  or  correction  of  any  errors  in  the  software  package  provided. 


6.  I/we  understand  that  no  material  subject  to  national  defense  security 
classification  or  proprietary  right  was  intended  to  be  released  to  us. 

I/we  will  report  promptly  the  discovery  of  any  material  with  such  restrictions 
to  the  Air  Force  approving  authority.  I/we  will  follow  all  instructions  con¬ 
cerning  the  use  or  return  of  such  material  in  accordance  with  regulations 
applying  to  classified  material,  and  will  make  no  further  study,  use  or  copy 
such  material  subject  to  security  or  proprietary  rights  marking. 


7.  I/we  understand  that  the  software  package  received  is  intended  for 
domestic  use  only.  It  will  not  be  made  available  to  foreign  governments 
nor  used  in  any  contract  with  a  foreign  government. 


Signature  of  Requestor 


Signature  of  Air  Force 
Approving  Authority 


Name/Title  of  Requestor  Name/Title  of  Air  Force 

Approving  Authority 


Organization/ Address 


Organization/Location 


City,  State  and  Zip  Code 
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LIST  OF  ACRONYMS  AND  SYMBOLS 


A 

<  A  > 

<  A'> 

A* 

AAM 

AB 

AFCMD/SA 

a 

C 

C3 

CHDF 

CSDF 

d 

erfc(x) 
ex,  ey 

F 

Gx>  Gy 

HOF 

h 

hI 

K 

<K> 

L 

M 
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Point  from  which  Penetrator  is  Attacked  ("Attack  Point") 
Expected  Attack  Point 

Expected  Intermediate  Attack  Point  (see  figure  11) 

Previous  Attack  Point 
Air-to-Air  Missile 
Air  Base  of  Interceptor 

Air  Force  Contract  Management  Division/Studies  and  Analyses  -- 
Office  Symbol  of  the  Directorate  of  Aerospace  Studies 

Characteristic  Distance  of  Extreme  Maneuver  (see  figure  5) 

Cosine  of  the  Total  Heading  (see  equation  (30d)) 

Command,  Control,  and  Communications 

Change-of-Heading  Distribution  Function 

Change-of-Speed  Distribution  Function 

Separation  between  Interceptor  and  Penetrator  (see  figure  8) 
Complementary  Error  Function 

Basic  Vectors  of  Air  Base-Centered  Cartesian  Coordinate  System 
(see  figure  7) 

Fence,  Tactical  Warning  (See  figures  1-3) 

Kernels  of  Convolution  Integrals  (see  equations  36)) 

Heading  Distribution  Function 
Penetrator  Cruising  Altitude 
Interceptor  Cruising  Altitude 
Point  at  which  Penetrator  is  Killed 
Expected  Point  of  Penetrator  Kill 
Expected  Penetrator  Path  (see  figure  3) 

Orthonormal  Transformation  Matrix  (see  equation  A7c)) 

Missile  Transit  Vector  (see  figure  12) 
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<Nk> 

NX>  Ny 

PDF 

P(t) 

Pit) 

<P(t)> 

PF 

Ph 

<P !  ( t )  > 

Pk 

<PK> 

Pt 

Q 

Q(x) 

q 

rbf 

rft 

ri 

Rl  ( t) 

rM 

RP 

rtr 


Expected  Number  of  Penetrators  Killed 
Normalization  Constants  (see  equations  (B15)) 

Position  Distribution  Function 

Penetrator  Position  at  Time  t 

Sensor-Observed  Penetrator  Position  at  Time  t 

Expected  Sensor-Observed  Penetrator  Position  at  Time  t 

Point  of  Fence  Penetration 

Location  of  Penetrator  Home  Base 

Expected  Interceptor  Position  at  Time  t 

Conventional  Kill  Probability 

Expected  Kill  Probability 

Location  of  Penetrator' s  Target 

Penetrator  Position  at  Start  of  Attack  (see  figure  7) 

Gaussian  Probability  Function 
Interceptor  Weapon  Load 

Distance  between  Penetrator' s  Home  Base  and  the  Fence  (see 
figure  1) 

Distance  between  Fence  and  Penetrator ‘s  Target  (see  figure  1) 
Total  Interceptor  Range 

Cumulative  Interceptor  Range  Expended  at  Time  t 
Operational  Range  of  Missile 
Total  Penetrator  Range 

Distance  between  Penetrator 's  Target  and  Penetrator 's 
Recovery  Base  (see  figure  1) 

Radius  Vector  to  Point  A 

Radius  Vector  to  Point  A* 

Raofus  Vector  to  Point  A’ 

,s  Vector  from  Air  Base  to  Penetrator  (see  figure  A2) 

•  i 


/ertcr  A:r  Base  to  Surveillance  Sensor  (see 


rc 

H) 

rSP 

s 

s 

SDF 

T 

-+■ 

Tl 

TSDF 

*A 

*A 

*F 

fcIN 

tK 

tT 

ta 

V[  v] 
V[v] 
V+[v+] 
v(t) 
v(t) 

A 

<v> 

V 

+ 

<v> 


Climb  Vector  (see  figure  9) 

Dive  Vector 

Radius  Vector  to  Point  Q  (see  figure  7) 

Radius  Vector  from  Surveillance  Sensor  to  Penetrator  (see  fi 
A2) 

Sensor 

Sine  of  Total  Heading  (see  equation  (30c)) 

Speed  Distribution  Function 
Penetrator  Transit  Vector  (see  figure  7) 

Interceptor  Transit  Vector  (see  figure  7) 

Total  Speed  Distribution  Function 

Attack  Time 

Projected  Attack  Time 

Time  at  which  Penetrator  Crosses  Fence 

Intermediate  Time  (see  equation  (97)) 

Kill  Time 

Time  of  Penetrator  Attack  against  US  Target 
Time  Interceptor  Arrives  at  Point  A 
Time  Interceptor  Leaves  Point  A* 

Total  Speed  Distribution  Function 
Speed  Distribution  Function 
Change-of-Speed  Distribution  Function 
Penetrator  Speed  at  Time  t 
Sensor-Observed  Penetrator  Speed  at  Time  t 
Expected  Sensor-Observed  Penetrator  Speed 
Maneuver-Induced  Change  in  Penetrator  Speed 
Expected  Maneuver-Induced  Change  in  Penetrator  Speed 


VC 

VD 


Interceptor  Climb  Speed 
Interceptor  Dive  Speed 
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vi 

vie 

Vm 

(*r»VV 

ivx,vy,vz) 

W 

WM 

WX>Wy 

xR(t)] 

^  + 

X-j  i X-j , X-j 

(x.y.z) 

(x.y.z) 
(<x>,<y>,<z  >) 

(xA,yA>zA) 

IxBP,yBp5ZBp) 

(xBS»yBS>2BS) 

(xp.yp.zp) 

(xq.yq.zQ) 

(xs,ys,zs) 

(xspiysp>zsp) 

(xT,yT,zj) 

Y[y(t)] 

Yi »Yi > Yi 


Penetrator  Speed  along  its  Expected  Path 
Interceptor  Speed 

Interceptor  Speed  along  its  Expected  Path 
Missile  Speed 

Penetrator  Velocity  Components  in  Spherical  Coordinates 
(see  equation  ( A7b ) ) 

Penetrator  Velocity  Components  in  Cartesian  Coordinates 
(see  equation  (A6b)) 

Auxiliary  Vector  in  Attack  Time  Solution  (see  equation  (57)) 

Auxiliary  Vector  in  AAM  Attack  Time  Solution  (see  equation 
(99b)) 

Cartesian  Components  of  W 
Position  Distribution  Function 

Terms  in  Convolution  Integral  Expansion  (see  equations  (37)) 

Cartesian  Coordinate  System  Centered  on  the  Interceptor  Air 
Base 

Sensor -Observed  Penetrator  Coordinates 
Expected  Sensor-Observed  Penetrator  Coordinates 
Coordinates  of  Attack  Point  A 
Coordinates  of  Attack  Point  A* 

Cartesian  Components  of  rgp 
Cartesian  Components  of  rgg 
Coordinates  of  Point  Pp 
Coordinates  of  Point  Q 

Cartesian  Coordinate  System  Centered  on  the  Surveil¬ 
lance  Sensor 

Cartesian  Components  of  rsp 

Coordinates  of  point  Pj 

Position  Distribution  Function 

Terms  in  Convolution  Integral  Expansion  (see  equa¬ 
tions  (37)} 


82 
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ac  ,a0 

8 

yc 

yd 

AX, Ay 
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Position  Distribution  Function 

Auxiliary  Terms  in  Flight  Profile  Calculations  (see 
equations  (75)) 

Angle  between  In  and  TI2  (see  figure  10) 

Interceptor  Climb  Angle 

Interceptor  Dive  Angle 

Auxiliary  Terms  in  Convolution  Integral  Expansion  (see 
equations  (37)) 

Auxiliary  Time  in  Attack  Time  Calculations  (see 
equation  (58)) 

Angular  Penetrator  Coordinates  with  respect  to  the 
Surveillance  Sensor  (see  figure  Al) 

Sensor-Observed  Angular  Penetrator  Coordinates 

Tangent  of  the  Sensor-Observed  Penetrator  Heading 

Tangent  of  the  Manuever-Induced  Penetrator  Heading 
Change 

Characteristic  Angle  of  Extreme  Manuever  (see  figure  5) 

Auxiliary  Term  in  Convolution  Integral  Expansion  (see 
equation  (32c)) 

Auxiliary  Terms  in  Position  Distribution  Functions  (see 
equations  (33)) 
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Standard  Deviations  in  Measurements  of  (r^e,^) 

Standard  Deviation  in  Penetrator  Speed  Measurement 

Maneuver-Induced  Uncertainty  in  Penetrator  Speed 

Standard  Deviations  in  Measurements  of  (Vr,  Vq.V^) 

Standard  Deviation  in  Total  Penetrator  Speed  (see 
equation  (26b)) 

Standard  Deviations  in  Measurements  of  (x,y,z) 

Auxiliary  Term  in  Position  Distribution  Function  (see 
equation  (36)) 

Standard  Deviations  in  Measurements  of  (x^p,  y$p,  z$p) 
Standard  Deviation  in  Penetrator  Heading  Measurement 
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Tc 

td 

Td 

tlag 

tl 

tm 

Ts 

tv 
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>p 
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Maneuver- Induced  Uncertainty  in  Penetrator  Heading 
Interceptor  Climb  Time 
Interceptor  Dive  Time 
Dead  Time 

Lag  Time  in  Interceptor 's  Transit  from  Point  A*  to  Point  A 
Induced  by  Flight  Profile 

Loiter  Time 

Missile  Flight  Time 

Search  Time 

Period  between  Successive  Updates  of  Vectoring  Information 

Heading  Distribution  Function 

Change-of -Heading  Distribution  Function 

Penetrator  Heading 

Sensor-Observed  Penetrator  Heading 

Expected  Sensor-Observed  Penetrator  Heading 

Maneuver- Induced  Change  of  Penetrator  Heading 

Expected  Maneuver-Induced  Change  of  Penetrator  Heading 
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